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The physics of low-energy quasi-particle excitations in disordered d-wave superconductors is a 
subject of ongoing intensive research. Over the last decade, a variety of conceptually and method- 
ologically different approaches to the problem have been developed. Unfortunately, many of these 
theories contradict each other, and the current literature displays a lack of consensus on even the 
most basic physical observables. Adopting a symmetry-oriented approach, the present paper at- 
tempts to identify the origin of the disagreement between various previous approaches, and to develop 
a coherent theoretical description of the different low-energy regimes realized in weakly disordered 
d-wave superconductors. We show that, depending on the presence or absence of time-reversal in- 
variance and the microscopic nature of the impurities, the system falls into one of four different 
symmetry classes. By employing a field-theoretical formalism, we derive effective descriptions of 
these universal regimes as descendants of a common parent field theory of Wess-Zumino-Novikov- 
Witten type. As well as describing the properties of each universal regime, we analyse a number of 
physically relevant crossover scenarios, and discuss reasons for the disagreement between previous 
results. We also touch upon other aspects of the phenomenology of the d-wave superconductor such 
as quasi-particle localization properties, the spin quantum Hall effect, and the quasi-particle physics 
of the disordered vortex lattice. 
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I. INTRODUCTION 

In recent years it has become clear that the super- 
conducting phase of the hole-doped cuprate supercon- 
ductors is of unconventional, d-wave, symmetry. Moti- 
vated by this observation, the properties of quasi-particle 
states in d-wave superconductors have come under in- 
tense scrutiny ]j]-pO|. A particular subject that turned 
out to be not only experimentally relevant but also theo- 
retically intricate is the influence of static disorder on the 
large-scale features of the d-wave phase. Leaving aside 
the pair-breaking effect of the disorder which strongly 
influences the high-temperature properties, attempts to 
resolve the impact of impurity scattering on low temper- 
ature spectral and transport properties of a disordered 
d-wave superconductor have ignited controversy in the 
recent literature. 

The key feature which distinguishes d-wave from the 
more conventional s-wave superconductors, and makes 
the disorder problem harder to solve, is the presence of 
four isolated "nodes" on the Fermi surface in the vicin- 
ity of which low-lying Dirac-type [^H quasi-particle exci- 
tations exist. Beginning with the work of Gorkov and 
Kalugin (lj, and later Lee 0, the earliest considera- 
tions, based on approximate self-consistent treatments, 
concluded that an arbitrarily weak impurity potential in- 
duces a finite density of states (DoS) at the Fermi surface, 
and leads to weak localization of all quasi-particle states 
in the two-dimensional (2d) system. Although the first 
result found support, at least superficially, in an analysis 
by Ziegler et al. Q, these conclusions were found to be 
in contradiction with those of other authors. 

In particular, following the early considerations of 
Refs. fH , Nersesyan, Tsvelik and Wenger § (NTW) 
developed a complementary approach in which the quasi- 
particle Green functions of a disordered d-wave supercon- 
ductor were mapped on a conformally invariant fermion- 
replica field theory. Based on this mapping, NTW pro- 
posed that the 2d system is characterized by critical prop- 
erties of the quasi-particle states at the Fermi level. Their 
model predicted that the DoS vanishes as v(E) ~ \E\ a 
with some disorder-dependent exponent a. 

To further add to this controversy, Senthil, Fisher, Ba- 
lents, and Nayak |(| (SFBN) proposed a description by an 



alternative fermion-replica field theory, namely the prin- 
cipal chiral non-linear sigma model over the symplectic 
group Sp(2r) (with r = 0). Void of any signatures of crit- 
icality, the large-scale behaviour of the two-dimensional 
model was argued to describe a phase of localized states, 
the "spin insulator" , thereby finding qualitative agree- 
ment with the early considerations and contradicting the 
findings of NTW. In contradiction to both NTW and the 
early works, the quasi-particle DoS at the Fermi level was 
found to vanish with a universal exponent of unity. 

Since this work, further analysis of the field-theoretic 
scheme generated fresh, and seemingly contradictory pro- 
posals. Fendley and Konik jl9j have argued that, when 
the scattering between pairs of opposite nodes is ne- 
glected, the low-energy properties of the weakly disor- 
dered system are described by a critical field theory that 
has Sp(2r) (with r — 0) for its target space and in- 
cludes a Wess-Zumino term. The same theory, now based 
on a peculiar implementation of a time-reversal sym- 
metry breaking lattice operator, had been suggested by 
Fukui ipHfl . Albeit belonging to a different universality 
class, both proposals pl| and [jl9| share key aspects with 
the theory of NTW, including the existence of extended 
quasi-particle states. 

How, if indeed at all, can the largely contradictory ap- 
proaches be reconciled with each other? And is it possible 
to unify elements of these approaches into one coherent 
theoretical formulation? This is the first complex of ques- 
tions we are going to address in this paper. 

A second and more pragmatic point we are going to 
address concerns the microscopic nature of the impurity 
scattering. It will turn out that much of the contro- 
versy in the field-theoretical literature described above 
can be attributed to the sensitivity of the d-wave quasi- 
particle physics to the range of the scattering potential. 
In a number of recent publications |^-^5|], it has been 
argued that both experimental and theoretical evidence 
hints at a strong enhancement of forward scattering in 
d-wave superconductors. Without going into details, we 
here merely state the essence of the argument, viz. that 
a significant amount of inter-node scattering would have 
a strong pair-breaking effect, thereby being in conflict 
with the very formation of stable d-wave order. Actually 
figuring out the meaning of the attribute "significant" 
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is a delicate issue we will not even attempt to address. 
Notwithstanding this uncertainty, one may ask the prin- 
cipal question of what kind of low-energy theory gov- 
erns the behaviour of quasi-particle Hamiltonians that 
are strongly anisotropic (albeit certainly not exclusively 
forward scattering) on the bare microscopic level. In 
particular, the question arises whether, as in disordered 
metallic systems, the impurity form factor merely leads 
to a renormalization of the scattering time or whether it 
might have a more substantial effect. As we are going to 
argue below, the latter is the case here: the extreme cases 
of pure forward scattering and isotropic scattering fall 
into distinct universality classes with qualitatively differ- 
ent properties. And, although the isotropic limit is ul- 
timately attractive (in the renormalization group sense) , 
one might speculate that under realistic conditions (fi- 
nite temperatures, experimental resolutions, etc.) the 
forward-scattering limit could well be the relevant one. 

A second aspect that makes the limit of weak forward 
scattering an interesting object of study concerns the 
physics of the mixed state: as with conventional type- 
II superconductors, the application of a magnetic field 
of intermediate strength (H c i < H < H C 2) forces vor- 
tices into the d-wave superconductor. The impact of vor- 
tex formation on the low-energy density of quasi-particle 
states of pure d-wave superconductors has been the sub- 
ject of an ongoing debate p6[-^l|. Recently, significant 
progress was made by Franz and Tesanovic [p8[ , who ex- 
ploited a novel singular gauge transformation, mapping 
the problem again on the problem of Dirac fermions. This 
transformation, which is analogous to that applied in the 
theory of the fractional quantum Hall effect, is applica- 
ble even in the presence of randomness. As a result, one 
obtains an effective low-energy model, essentially Dirac 
fermions with random scalar and vector potential, that 
can be applied to explore the quasi-particle DoS of the 
disordered vortex lattice. We will return to this subject 
below. 

Although, for a non-expert, the detailed field- 
theoretical analysis contained in this paper is often some- 
what involved, we believe that the main conclusions of 
our survey of the dirty "d-wave system" — for lack of 
a better terminology we will often refer to the quasi- 
particles of a disordered <i-wave superconductor by that 
word — should be widely accessible. We have, therefore, 
chosen to summarize the complete phase diagram here 
in the introduction. At the same time, this allows us 
to organize and place into context many of the existing 
field-theoretical works on the subject. 



A. Intrinsic and Extrinsic Symmetries 

While, at first sight, the various points brought up 
above seem unrelated, the microscopic analysis below will 
show that they do, in fact, all find their origin in a com- 
mon microscopic mechanism. Just as for the clean d-wave 



Gorkov Hamiltonian |3^,|33|], the Hilbert space for the 
low-energy Hamiltonian of the pure forward-scattering 
system foliates into four sectors not connected by the 
disorder potential. These nodal sectors can be grouped 
into pairs related by parity. The low-energy physics of 
the individual sectors turns out to be more intricate than 
that of the complete system, which is comprised of four 
Dirac nodes coupled by impurity scattering. In general 
it will, of course, be the behaviour of the full system that 
matters. There are, however, a number of situations in 
which signatures of the decoupled system remain visible. 
Exploring these scenarios, which are accompanied by a 
drastic change in the phenomenology of the system, is 
one of the major issues to be addressed in this paper. 

The origin of the relative non-triviality of the decou- 
pled system can be understood qualitatively by noticing 
that the structure of any theory describing the low-lying 
quasi-particle excitations of a disordered system is es- 
sentially fixed by two types of symmetries: "intrinsic" 
symmetries, such as the behaviour of the Hamiltonian 
under spin-rotation, time-reversal, particle-hole transfor- 
mations etc., and "extrinsic" symmetries such as transla- 
tional or rotational invariance or, for that matter, parity. 
Specifically, the unperturbed d-wave superconductor is a 
particle-hole symmetric, spin-rotation and time-reversal 
invariant system, thereby falling into the symmetry class 
CI in the classification of Ref . p4fl . The implementation 
of these intrinsic symmetries into an effective low-energy 
theory of a superconductor system was formulated some 
time ago by Oppermann |35| l , and was recently rederived 
in the context of d-wave superconductivity by SFBN || . 

The behaviour of the d-wave system under extrinsic 
symmetry operations is more remarkable. In fact, as we 
will see, it is disregard for these symmetries that has 
driven astray some of the field theories based on intrinsic 
symmetry: not all terms allowed by the intrinsic sym- 
metries respect the extrinsic symmetries. The point is 
that individual nodes are non-invariant under "parity" , 
by which we mean any operation that reverses the orien- 
tation of two-dimensional space |56|. Of course, the full 
four-node system emerging from a manifestly invariant 
Gorkov Hamiltonian is parity invariant, but individual 
nodes are not. Remarkably, this parity non-invariance 
and the intrinsic symmetries conspire to let the system 
of isolated nodes fall into a symmetry class that differs 
from the Wigner-Dyson classes and the classes commonly 
attributed to bulk superconductors. The complete phase 
diagram of the d-wave system in fact separates into a to- 
tal of four distinct regions (see Fig. |l|) distinguished by 
the presence/ absence of time-reversal invariance and by 
the correlation radius of the disorder potential. Taking 
each region in turn, we here summarize the phenomenol- 
ogy and the basic properties of the respective low-energy 
theories jjTj. 
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Hard Scattering 



CI Spin Insulator 



AIII Critical Metal 



C (SQHEI) 

► T-Breaking 

A (SQHEII) 



FIG. 1. Symmetry classes realized in dirty d-wave super- 
conductors with conserved spin. The two types of perturba- 
tion causing crossover between classes are indicated. 



B. Symmetry Class AIII 



So, beginning with the time-reversal invariant system 
with only forward scattering, one finds that the corre- 
sponding symmetry class is not CI, but rather AIII, i.e. a 
symmetry class that is typically realized in disordered rel- 
ativistic (chiral) fermion systems. Indeed, it is straight- 
forward to verify that the low-energy physics of individ- 
ual nodes is described by a model of Dirac fermions in 
the presence of some quenched random gauge field. The 
implications of this correspondence for the disordered 
system were first noticed and analysed in the seminal 
work of NTW [||, who pointed out that the low-energy 
physics of the individual nodes is critical. NTW argued 
that a general mechanism uncovered long ago by Wit- 
ten ]3j| applies in particular to the <i-wave system: in a 
Lagrangian formulation based on non-Abelian bosoniza- 
tion, non-trivial transformation behaviour with respect 
to parity is reflected by the presence of a Wess-Zumino- 
Novikov-Witten (WZW) term in the action functional. 
Specifically, NTW mapped the d-wave system on a field 
theory with U(l) x SU(r) x SU(iV) symmetry (r being 
the number of replicas). In the limit of isolated nodes 
(N = 1), the theory has the field manifold U(l) x SU(r), 
with the SU(r) sector being governed by a so-called level- 
1 WZW action. The latter theory enjoys the existence 
of an infinite-dimensional (current algebra) symmetry, 
which makes it solvable to a large extent. Several of its 
features, including the existence of algebraically decaying 
correlations (absence of localization!), and the scaling of 
the DoS in the vicinity of the Fermi energy can be derived 
rigorously. 

Within the framework of the supersymmetric formu- 
lation, the properties of the forward-scattering system 
are found to be described by two copies of a WZW 
model of level 1 and type A\A, i.e. a WZW model 
whose fields take values in a supermanifold that, for rea- 
sons discussed in body of the paper, is labeled by A\A. 
Each copy derives from a pair of nodes on the Fermi 
surface related by parity or inversion symmetry. For- 
mally, the total zero-energy effective action is given by 
S = W[Mi,7] + WlMi 1 ,*/- 1 ] + g' Jd 2 rO' 00 where the 
dimensionlcss parameter 7 (defined microscopically later) 
specifies the degree of anisotropy of the Dirac cones, and 



W[M,i\ = JLr[M\+g J d 2 rO 0Q 

— [ d 2 rSTr {^dxM^dxM + -id 2 M- 1 d 2 M) 

87T J 



is the anisotropic WZW action. The fields Mi and M 2 
take values in (a maximal Riemannian subspace of) the 
supergroup GL(2|2), Y[M] represents the WZW term, 
and Ooo, Cqo are low-angle scattering operators whose 



specific structure will be described in Section IV. Note 
that, since the WZW term is odd under M — ► M -1 
(i.e. T[M- X ] = -T[M]), the WZW terms for the two 
sectors of nodes carry opposite sign, which reflects the 
fact that they derive from Dirac operators carrying op- 
posite orientations. Further, the fact that GL(2|2) and 
not GL(1|1) - the supersymmetric analog of a fermion 
replica theory with group manifold U(r) - appears as 
the relevant degree of freedom, has to do with the be- 
haviour of the system under time reversal. This aspect, 
unimportant for the physics of the system of isolated 
nodes, becomes vitally relevant once inter- node scatter- 
ing is switched on. 

Perhaps the most characteristic feature of models of 
Dirac fermions subject to Abelian gauge randomness is 
their stability under renormalization. Indeed, the low- 
energy theory defined by the action W is a fixed point 
of the renormalization group for each value of the di- 
mensionless disorder strength g, as long as g does not 
exceed some critical value. The absence of running cou- 
plings bears consequences for various physical observ- 
ables. One finds that, for low energies, the DoS vanishes 
as v(E) oc \E\ a , with the g-dependent family of expo- 
nents 

1 - 2g A 

a = — , 

1 + 2 9 /tt 

in accord with the analysis of NTW. We also recapitulate 
that the eigenstates of the random gauge Dirac system 
are extended [j39[ . 

While exhibiting critical spectral and transport prop- 
erties, the effective action for c£-wave superconductors in 
class AIII is unstable against the influence of short-range 
isotropic scattering, a matter not correctly handled by 
NTW. Impurity scattering involving large momentum 
transfer couples the nodes and leads to a "locking" of 
the fields at the different nodes. This process bears two 
consequences for the structure of the low-energy theory: 
first, and in contradiction to the naive expectation, the 
field configurations surviving the locking of the formerly 
independent GL(2|2) fields take values in a different man- 
ifold, viz. OSp(2|2). Second, the locking Mi = M 2 im- 
plies that in the full action the WZW terms cancel, which 
makes the restoration of parity invariance in the coupled 
system explicit. (In parentheses we note that the absence 
of a WZW term is a sufficient but not a necessary con- 
dition for parity invariance of physical observables: for a 
single WZW theory, the sign change in T[M] caused by 
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a parity operation can be compensated by a field trans- 
formation M = e x — > M^ 1 = e~ x , i.e. the field X is 
a pseudo-scalar. However, in a system comprised of two 
coupled sub-theories carrying opposite orientations, this 
compensating transformation is inhibited on cither sub- 
theory. Moreover, there exists no room for a WZW term 
in the locked full theory, as the sum of coupled Dirac 
theories for the d-wave system is anomaly- free.) Since it 
was the WZW term that rendered the NTW theory crit- 
ical, it is evident that the low-energy behaviour of the 
full theory will be of different type. 

C. Symmetry Class CI 

A fermion-replica theory of the full system, a non- 
linear sigma model on the field manifold Sp(2r), was first 
derived and analysed by SFBN. Essentially, this model is 
a ci-wave version of the Sp(2r) non-linear sigma model 
formulation of time-reversal invariant superconducting 
glasses developed earlier by Oppermann |35[ ). Void of 
any elements supporting critical behaviour, the predic- 
tions of that theory differ greatly from NTW: all quasi- 
particle states are localized, a fact that led SFBN to coin 
the term "spin insulator" ; the DoS still vanishes at zero 
energy but in a manner altogether different from the scal- 
ing obtained by NTW. All in all, the low-energy theory 
of the composite four-node system seems to have little 
overlap with the critical theory of the decoupled forward- 
scattering system. 

Once the fields belonging to the different nodal sectors 
become locked, the total effective action from above takes 
the form S[M] = W[M^] + W[M _1 ,7 -1 ], i.e. 

S[M] = - 7 +J 1 J d 2 r STr (fyM-^A/) , 

where the field M now belongs to the group manifold 
OSp(2|2). Technically this is a principal chiral non-linear 
sigma model of type D\C '. It represents a supersymmetric 
extension of the action derived by SFBN M , and leads to 
the prediction of quasi-particle localization of all states 
in two dimensions — the spin insulator phase — and a 
DoS that vanishes at zero energy. 

It seems unlikely that symmetry class Alll is realized 
in the physical d-wave system. As mentioned earlier, the 
bare microscopic theory will probably contain a large 
forward-scattering amplitude, plus a certain amount of 
inter-node scattering rendering it a member of symme- 
try class CI. This hard-scattering admixture is a relevant 
perturbation which means that the theory will ultimately 
flow towards the strong-coupling fixed point in class CI. 
On the other hand, realistic experimental conditions are 
likely to prevent the system from exploring all limiting re- 
gions of the phase diagram. This motivates us to ask how 
fast the flow towards the strong-coupling regime takes 
place and what properties the crossover theories, situ- 
ated somewhere in between the independent-node and 



the strongly coupled theory possess. We will return to 
this issue below. 



D. Symmetry Class C 

An alternative mechanism of allowing elements of the 
decoupled theory to resurface is by breaking parity ex- 
plicitly. At least theoretically this is realized in the so- 
called spin quantum-Hall (SQH) system [ ^2pfj[| , where 
the nodal degeneracy is lifted via the addition of an id xy 
component to the order parameter. The resulting low- 
energy phase shares much of its phenomenology — the 
formation of edge states, quantized values of the trans- 
verse component of the conductivity tensor, etc. — with 
the integer quantum Hall effect and has, therefore, been 
christened the spin quantum Hall state. (The prefix 
"spin" indicates that in the superconductor problem, the 
behaviour of spin currents or thermal currents is consid- 
ered. Because of the non-conservation of quasi-particle 
charge, the study of charge transport coefficients becomes 
meaningless.) Originally, it was suggested j4(]] that the 
formation of a secondary id-order parameter component 
giving the d + id state might be quasi-spontaneously 
driven by a weak external magnetic field. Recent work 
[fflf has cast some doubt on the applicability of this sce- 
nario and, therefore, on the experimental relevance of 
the spin quantum Hall effect. Theoretically, however, the 
d+id system remains an interesting object of study: first, 
it represents a rare example of a system with an exactly 
solvable quantum-Hall phase transition driven by disor- 
der E4). Second, unlike with d-wave superconductors in 
the previously discussed symmetry class CI, aspects of 
the physics of individual nodes survive the thermody- 
namic limit. 

In fact, the id-component does more than lift the node 
degeneracy, implying that it does not lead to a resurrec- 
tion of the WZW model: the inclusion of an id xy compo- 
nent not only breaks parity but also time-reversal invari- 
ance, and hence reduces the symmetry of the system from 
CI down to C. Within a fermion-replica approach, the 
class C system is described by a field theory taking val- 
ues on the manifold Sp(2r)/U(r). This field theory does 
not support a WZW term. It does, however, admit for 
the existence of a closely related parity-breaking term in 
the Lagrangian, viz. a topological theta term. Based on 
phenomenological considerations, it has in fact been ar- 
gued that a non-linear sigma model over Sp(2r)/U(r) 
with a theta term should be the relevant theory of the 
disordered d + id system. 

Exploring how such a theory may emerge as a descen- 
dant of the class CI theory of the non-degenerate node 
system will be a further topic to be addressed below. The 
explicit form of the resulting supersymmetric soft-mode 
action for the system of class C reads as 

-| J d 2 r STy (a^Qd^Q + <j° 12 e^Qd^Qd,Q) , (1) 



5 



where the WZW term has transformed into a topologi- 
cal theta term. In contrast to the conventional integer 
quantum Hall effect exhibited by normal systems, the 
fields of the saddle-point manifold Q belong to the coset 
space OSp(2|2)/GL(l|l), with the relevant Riemannian 
symmetric superspace inside that manifold being of type 
Dili | CI. 

Unlike its normal relative, the critical theory for the 
spin quantum Hall effect for class C is readily amenable 
to theoretical investigation. Indeed the properties of this 
phase have been explored in a number of recent numeri- 
cal and analytical works. Surprisingly, a mapping pl| of 
the class C theory onto a network model |^3j shows the 
transition to be in the same universality class as classical 
percolation jl2| . This correspondence, which applies only 
for the statistics of various low-order moments, relies on 
a remarkable cancellation of interference channels, and 
allows several of the critical exponents to be determined 
exactly p^Ji^| . Very recently, a proposal for the gen- 
eral critical theory of the spin quantum Hall transition 
has been made by Bernard and LeClair , although its 
validity still remains a matter of some controversy. 



E. Symmetry Class A 

To exhaust the number of different low-energy scenar- 
ios realized in the cZ-wave system, one may consider a sit- 
uation where a perturbation breaking time-reversal sym- 
metry is superimposed on an impurity potential causing 
only forward scattering. Below we will identify the re- 
sulting system as a member of symmetry class A, i.e. 
the standard Wigner-Dyson class of unitary symmetry. 
We find that, as with its class C counterpart, the class 
A system also supports a (spin) quantum Hall transi- 
tion. Somewhat surprisingly, however, that transition 



turns out to fall into the universality class of the stan- 
dard integer quantum Hall transition. 

Formally, the low-energy theory is again described by 
an action functional including a theta term, similar to 
the one above. The main differences with the previ- 
ously discussed case are that (a) one obtains two copies 
S[Q,j] of the action, one for each nodal sector, and (b) 
that the fields Qi (i = 1,2) belong to the coset space 
GL(2|2)/(GL(1|1) x GL(1|1)). This means that the field 
theory is of type AHI|.AIII which is the supersymmetric 
extension of Pruisken's theory for the conventional inte- 
ger quantum Hall effect. 

Finally, we remark on a second and, arguably, more 
physically motivated scenario in which the classes A and 
C are engaged — the mixed phase of a type-II supercon- 
ductor. On symmetry grounds alone, one would expect 
the Hamiltonian of the vortex phase of the dirty d-wave 
system to belong to one of these classes. Although the 
nature of the quasi-particle states in the vortex phase of 
the clean cZ-wave system has been addressed in a num- 
ber of (seemingly controversial) publications 1 2G 3l]|I(| . a 
field-theoretic analysis of this disordered system has not 
been given in the literature. This motivates us to focus 
our discussion of the symmetry classes C and A on the 
mixed state, not on the spin quantum Hall system. 

Furthermore, for brevity, we have chosen not to com- 
ment upon the influence of the breaking of spin rotational 
invariance either by magnetic impurities or spin-orbit 
scatterers. Such systems, which lie outside the scope of 
present work, have been addressed in the recent litera- 
ture and we refer to Refs. |47]-|4^] for a discussion of the 
phenomenology of these systems. 

Altogether, the key characteristics of the field theories 
for the d-wave system are summarized in the following 
table: 



Class 


T 


hard scattering 


Goldstonc modes 


saddle-point manifold 


NLctM 


fcrmion-rcplica analog 


Am 


+ 




GL(2|2) x GL(2|2) 


GL(2|2) 


A\A 


U(2r) 


A 






GL(22) 


GL(2 2)/(GL(l|l) x GL(1|1)) 


,4111 1 ,4111 


U(2r)/(U(r) x U(r)) 


CI 


+ 


+ 


OSp(2|2) x OSp(2|2) 


OS P (2|2) 


D\C 


Sp(2r) 


c 




+ 


OSp(2 2) 


OS P (2 2)/GL(l|l) 


L>III|CI 


Sp(2r)/U(r) 



Note that the contents of this table are in complete agree- 
ment with Ref. ]50| ], where each of the ten symmetry 
classes, as enumerated by the large families of Cartan's 
list of irreducible symmetric spaces, was put in corre- 
spondence with the saddle-point manifold of a non-linear 
sigma model. This coincidence is by no means acciden- 
tal. Although Ref. (5Cj was formulated in the random- 
matrix setting, the translation scheme employed there is 
of universal validity and unerringly identifies the internal 
(field-theoretical) global symmetries for each symmetry 
class. The present treatment has the virtue of being well 
adapted to the specific case of d-wave superconductors, 



and shows quite explicitly how the hierarchy of Goldstonc 
modes emerges in this context. 



F. Self-consistent Theories 

Before leaving this introductory section, let us re- 
turn to the controversy between the predictions made 
by the various field-theoretical approaches, and the self- 
consistent diagrammatic schemes developed in the liter- 
ature. Readers who find the following, purely qualita- 
tive, remarks cryptic are referred to Section M for a more 



G 



substantial discussion. The last few years have seen a 
build-up of confusion concerning the behaviour of the 
density of states at very small energies: straightforward 
diagrammatic analysis predicts a finite DoS at zero en- 
ergy, whereas, in the field theory approaches, unbounded 
fluctuations of Goldstone modes lead to a vanishing DoS. 
An intuitive argument in favor of a finite DoS is that any 
average over random potential fluctuations should intro- 
duce some imaginary self energy £ in the Green function 
of the Dirac operator. This self energy would lead to a 
smearing of spectral structures over some energy window 
of width given by £ and would, therefore, be in conflict 
with the formation of a cusp in the energy-dependent 
DoS. This argument, however, is overly simple; in par- 
ticular it does not account for the full complexity of the 
scattering problem in random Dirac (or gapless super- 
conductor) systems. Let us try to discuss this point on a 
somewhat more technical level: 

The various classes of diagrams appearing in pertur- 
bative approaches can be grouped, roughly speaking, 
into two categories: diagrams of the self-energy type 
(commonly evaluated in SCBA or similar approxima- 
tions), and various types of ladder diagrams (alias dif- 
fusion modes). While the former are usually associ- 
ated with some smearing of average single-particle quan- 
tities by disorder, the latter describe mesoscopic fluc- 
tuations which become operative at long range. More 
precisely, ladder diagrams describe the long-range quan- 
tum interference between retarded and advanced single- 
particle Green functions. For normal-conducting sys- 
tems, the distinction between these diagrammatic ele- 
ments is perfectly canonical: single-particle Green func- 
tions are renormalized by self-energy diagrams whereas 
higher-order Green functions with poles on both sides 
of the real energy axis may be influenced by diffusion 
modes. 

For superconductors however, and for gapless su- 
perconductors in particular, this distinction breaks 
down. In previous diagrammatic approaches, the DoS of 
dirty d-w&ve superconductors was computed with (self- 
consistent) account for self-energy diagrams. This led to 
some smearing and, therefore, to a finite DoS. However, 
as was shown in Ref. H] for the prototypical case of su- 
perconductor/normal quantum dots, the single-particle 
DoS of a gapless superconductor may be affected by 
diffusion-like modes. This phenomenon can be under- 
stood on various levels of sophistication. It can be ex- 
plained by qualitative semi-classical reasoning [B4| or, 
more technically, by the fact that unlike in normal sys- 
tems the single-particle Green function of a supercon- 
ductor possesses poles on both sides of the real energy 
axis. While for a bulk s-wave superconductor, all poles 
are shifted far into the complex plane by the bulk or- 
der parameter (a manifestation of Anderson's Theorem 
pl[), the situation in gapless systems is different. Here, 
multiple impurity scattering suppresses the self energy. 
Technically, the suppression is caused by ladder diagrams 
that sum up to yield diffusion modes. These modes enter 



already at the single-particle level |34j] , and are missing 
from previous diagrammatic works on d-wave supercon- 
ductivity. 

In the limit of zero energy, the diffusion modes become 
massless (or of infinite range) implying that a purely 
perturbative approach is met with a problem. (Tech- 
nically, in two dimensions, the contribution of each diffu- 
sion mode scales as g^ 1 \u(E/Eq), where E$ is some cut- 
off energy and g a measure of the dimensionless (spin) 
conductance of the system. For E sufficiently small, 
the expansion becomes uncontrolled.) At this point it 
is good to remember that all the concepts discussed 
above have a ficld-thcoretical analog, with an optional 
non-perturbative extension. In particular, the diffusion 
modes of the diagrammatic approach have the signifi- 
cance of perturbative excitations of the Goldstone modes 
discussed above — very much like a spin wave can be 
interpreted as a perturbative manifestation of the Gold- 
stone mode related to the intrinsic spin-rotation symme- 
try of a ferromagnet. The fluctuations of these modes, 
restricted only by finite energies, become stronger and 
stronger as the energy approaches zero. For sufficiently 
small energies, (a) perturbative (diagrammatic) schemes 
for treating these modes are ruled out, and (b) the non- 
perturbative analyses discussed below produce a vanish- 
ing DoS. (In passing we note that such mechanisms are 
not unfamiliar in mesoscopic physics. For example, for 
energy differences u> — E\ — E^ larger than the mean level 
spacing A, spectral correlations of generic normal disor- 
dered systems can be treated by diffusion-mode diagram- 
matic approaches. For to/ A < 1 perturbation theory 
breaks down, a phenomenon that manifests itself in the 
appearance of unphysical divergences. In this regime, a 
non-perturbative integration over the appropriate Gold- 
stone mode |52| produces the correct result.) Summa- 
rizing, we find that (i) the issue of the DoS is crucially 
related to long-range modes of quantum interference op- 
erating on the single-particle level, (ii) to obtain the cor- 
rect, vanishing, behaviour of the DoS, a non-perturbative 
treatment of these modes in inevitable but that (iii) sig- 
natures for the incompleteness of previous perturbative 
analyses can be found within the diagrammatic frame- 
work. 

Now, it is important to realize that the arguments 
above by no means amount to a categoric prediction of a 
vanishing zero-energy DoS in bulk d-wave superconduc- 
tors. The point is that the weak disorder field-theoretical 
approach is based on certain model assumptions which 
leave room for complementary scenarios. Therefore, to 
complete our introductory discussion, let us briefly men- 
tion a number of competing theories of the disordered d- 
wave system which lie outside the field-theoretic scheme 
outlined above. Typically such theories rely on pecu- 
liarities of the d-w&ve system. For example, Balatsky 
and Salkola j^] have proposed a mechanism whereby the 
weak coupling of "marginally-bound" impurity states at 
the Fermi level leads to the absence of quasi-particle lo- 
calization. 
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When subject to a single impurity in the unitarity 
limit, a quasi-particle state is created at zero energy 
|E^,E3. In a recent work, Pepin and Lee (l^] have pro- 
posed that, when subjected to many such impurities, 
these quasi-particle states broaden into a narrow delo- 
calized band around zero energy. On this basis, these 
authors have proposed that the properties of the dirty 
d-wave system at very low energies are characterized by 
extended quasi-particle states. Based on the RG analy- 
sis contained in this paper, it is our belief that, even for 
scatterers in the unitarity limit, there must ultimately be 
a cross-over at low-energy scales to the physics advocated 
by SFBN. 

Finally, some authors pcfl have inquired into the im- 
portance of maintaining self-consistency of the d-wave 
order parameter in a disordered background. As others 
before, in this work we will take a pragmatic view and 
will not attempt a synthesis of self-consistent aspects of 
the theory. 

A second and probably more severe limitation regards 
the role of interactions. As with more conventional dis- 
ordered normal and superconductor systems, the low- 
energy physics of quasi-particles in d-wave supercon- 
ductors, too, will be affected by various mechanisms of 
inter-particle interactions: zero-bias anomalies, interac- 
tion contributions to transport coefficients, interaction 
induced renormalization of the propagation in the Cooper 
channel, and more. Beyond that, one should expect that 
collective excitations, both of purely electronic and of 
phononic type, contribute to (thermal) transport coeffi- 
cients, and further modify the low-energy behaviour of 
quasi-particle states. None of these effects will be con- 
sidered in the present paper. (For a discussion of inter- 
action effects in d-wave superconductors, see Ref. [[l4j|.) 
Needless to say, this may impose severe limitations on 
the relevance of the present work for the physics of real 
d-wave superconductors. 

To what extent do the predictions of the weakly disor- 
dered theory bear comparison with experiment? As men- 
tioned above, superconductors of unconventional d-wave 
symmetry are realized in the high-temperature cuprate 
superconductors. Here the influence of both magnetic 
and non-magnetic impurities has been investigated in- 
tensively. In d-wave systems, both types of disorder are 
pair-breaking and lead to a rapid depression of T c with 
doping. However, at very low doping concentrations, T c 
remains a significant fraction of its optimal value. In this 
case the influence of disorder on the low-energy quasi- 
particle properties have been studied extensively. 

In particular, recent angle resolved photoemission ex- 
periments |55| confirm both the integrity of the large- 
scale nodal structure of the spectrum and confirm the ex- 
istence of long-lived quasi-particle states near the Fermi 
level. The same measurements show a large anisotropy 
of the Dirac nodes with, for example, Bi2Sr2CaCu20s 
exhibiting a ratio of t/A = 7 = 20. 

At very low temperatures, phase coherence properties 
of the quasi-particle states should be accessible via spin 



transport. The latter have been investigated indirectly 
through measurements of the thermal conductivity at 
very low temperatures: relating the thermal conductance 
k(T) to the quasi-particle spin conductivity a s through 
the Einstein relation, one obtains the Wiedemann-Franz 
ratio k(T)/To s = (4-7r 2 /3)fc|. Taking the bare value for 
the spin conductance from theory, one obtains p,[15[ 



k(T«1/t) _n 



= 4—ki 



Remarkably, this result is found to be in close agreement 
with experiment |57j . In BiSCCO the extrapolated value 
of the longitudinal part of the thermal conductivity re- 
mains constant over a wide range of impurity concentra- 
tions. This result has been triumphed (see, e.g., Ref. [[38)) 
as experimental support for the self-consistent theory of 
Lee §. 

However, it should be noted that this result does not 
sit comfortably with the class CI theory, which predicts 
that the spin conductance should exhibit weak local- 
ization corrections which renormalize a s down from its 
universal bare value. Whether the seeming coincidence 
of the measured conductance and the bare value signi- 
fies strongly enhanced forward scattering, or whether, 
as seems more likely, it can be explained in terms of 
a short phase coherence length, is as yet unclear. On 
this point, it is interesting to note that low-temperature 
measurements of the quasi-particle conductivity in the 
underdoped cuprate YBa2Cu40s |5^] indicate that the 
quasi-particle states are localized. Although it would be 
useful to explore the thermal Hall transport |30| , exper- 
iments have been so far unable to access the same very 
low temperature regime. 

The paper is organized as follows: in Section [n] we will 
formulate some basic concepts required for the construc- 
tion of a low-energy theory of the d-wave system. This 
will include a symmetry analysis of the disordered Hamil- 
tonian, its representation in terms of a (supersymmetric) 
generating functional, and a renormalization group anal- 
ysis of the relevancy of the various channels of disorder. 
In Section III we ask how an effective low-energy theory 



can be distilled from the microscopic functional integral. 
We will find that standard schemes, such as a straightfor- 
ward gradient expansion, are met with severe difficulties. 
In Section IV, non-Abelian bosonization will be applied 



as an alternative method to derive the low-energy theory. 
In a number of subsections the predictions of this theory 
for the four symmetry variants introduced above will be 
discussed. In Section we compare these results with 
the phenomenology predicted by other approaches, such 
as self-consistent T-matrix approaches. To complement 
the theoretical analysis of the low-energy properties of 
the d-wave system, in Section [v| we present a limited re- 
view of the progress made in the numerical investigation 
of the density of states and show how these results relate 
to the different regimes discussed in the tex t. Finally, we 
close with a brief discussion in Section VII. 
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II. DIRTY D-WAVE (BASICS) 

In this section we discuss the essential concepts needed 
for the construction of a low-energy effective theory de- 
scribing the quasi-particles of a disordered spin-singlet 
superconductor with d x 2_ y 2 symmetry. We start out by 
reviewing some elementary features of the Gorkov Hamil- 
tonian for that system. Much attention will be paid to 
the generic symmetries, which turn out to be quite sub- 
tle. We then express the Green functions of the Hamil- 
tonian as a field integral. The result is a supersymmetric 
generating functional, which will serve as the basic plat- 
form for the low-energy effective field theories developed 



in later sections. Similar considerations for a class of re- 
lated two-dimensional systems, having a low-energy limit 
modeled by a disordered Dirac equation, were first made 
by Fradkin [p2[ . 

A. The Hamiltonian 

Basic to the theory of quasi-particle transport in spin- 
singlet superconductors are the Gorkov equations [p2| . 
They follow from a quasi-particle Hamiltonian which in 
second quantized language, and on a lattice of sites i, is 
written as 




where are the hopping matrix elements, and Ay is the 
lattice d-wave order parameter. Time-reversal and spin- 
rotation invariance constrain these matrices to be real 
symmetric: tij = tji = iy and Ay = Aji = Ay. Their 
precise form will be specified shortly. The operators c\ a 
create spin- 1/2 fermions, /i is the chemical potential, and 
the sum 53 (y) ex tends over nearest neighbours on a two- 
dimensional square lattice with spacing a. By making a 
particle-hole transformation on down spins, 

d T = c T , d l = c\ , 

the Hamiltonian is cast in the concise form 

H =} d}((tij - ^Sij)a 3 + Aijai)dj , 

where a\, 02, 03 are the Pauli matrices. By way of warn- 
ing, we remark that this transformation makes it seem 
as though charge were conserved and spin was not. In 
reality, the opposite is true: quasi-particle charge is not 
conserved, whereas spin is. To facilitate future manipu- 
lations, we perform a 7r/2-rotation d 1— > exp(i7rcri/4)d, to 
transform to 

H = J2 d l( (^j ~ )02 + A tf o-i )dj . (3) 

(ij) 

The matrix elements iy and Ay consist of clean and 
dirty parts. The clean part of the Hamiltonian, Hq, 
which we focus on first, assumes its simplest form in mo- 
mentum space: 

Ho=J2 d l [(*(*) - M>a + A(*)<ri] d k . 

k 

The kinetic energy and the d-wave order parameter are 
taken to be 



-t^+yhi ) ( ct] ) ' (2) 

1 

t{k) = t(cos(k x a) + cos(k y a)), 
A(k) — A(cos(k x a) — cos(k y a)). 

It is the momentum dependence of the order parame- 
ter A(k) that distinguishes c?-wave superconductors from 
their relatives with s-wave pairing. The dispersion rela- 
tion of H , 

e(k) = ±VWfc)-^) 2 + A(fc) 2 , 

is displayed in Fig. || for the special value fj, = (the half 
filled band). 




FIG. 2. The quasi-particle spectrum of a clean d-wave su- 
perconductor showing the existence of four Dirac nodes. 

The most important feature of the function e(k) is that 
it vanishes (for fi — 0) at the four points 

(k x ,k y ) = a - 1 (±7r/2,± 7 r/2) . 

When the chemical potential is shifted away from half 
filling, these nodal points move to different locations in 
the Brillouin zone. No qualitative change is caused by 
that, and we shall therefore stick to the case /i = 0. Af- 
ter the introduction of small momentum offsets by setting 
(k x ,k y ) — (±7r/(2a) + 5k x , ±n/(2a) + Sk y ), the disper- 
sion relation in the vicinity of the nodal points is approx- 
imated as 
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e(k + 8k) ~ ±a 



t 2 H 



A 2 k 2 , 



where ki = 5k x =p 5k y . Thus the dispersion relation near 
the nodal points is linear in k, and the low-energy physics 
will therefore be of Dirac type. Anticipating that we will 
be primarily interested in the nodal regions, we introduce 
four species of fermion operators, 



k n +8k , 



where n G {1,1,2,2} labels the nodes, with the assign- 
ments to node momenta k n being 

!_> ^1(^/2,^/2), 
2 -> a" 1 (7r/2,7r/2) , 
S^a-^-TrA-Tr^) , 
1 -> a _1 (7r/2,-7r/2) . 



The overbar over 8^ does not denote complex conjuga- 
tion, but is motivated by the observation that, for 7=1, 
vanishes on holomorphic functions f(x% +1x2), while 
c^ 1 -* annihilates antiholomorphic functions f{x\ — ix 2 ). 

We now introduce disorder. In general, both the nor- 
mal part of the Hamiltonian and the order parameter will 
contain a random component, which we denote by Vi mp 
and Aj mp , respectively. Specifically, we assume that 



Vi 



imp 



and A 



imp 



are uncorrelatcd Gaussian dis- 



tributed variables with 

• zero mean: (V imp ) = (A imp ) = 0, 

• and the same variance: 



(X imp (r)X imp (r')) = |/(|r-r'|/0: 



If we restrict attention to the low-energy sector, the 
Hamiltonian of the clean system decomposes into four 
nodal sub-Hamiltonians: 



k ■ 



n , k 



With the definition of the characteristic speed v and the 
anisotropy parameter 7 by 



ta 



t 

1= A' 



the continuum real space representations of the operators 



Hg are 



Hi = -Hi = -w( 7 - 1 <7 1 9 1 + a 2 d 2 ) 



H, 



-Hq = —iv(y 1 a 1 d 2 + cr 2 di) . 



Here di = d Xi [i = 1, 2) and xi = x =p y. 

The structure of Hq makes the connection between 
clean d-wave superconductors and the two-dimensional 
Dirac Hamiltonian manifest. An equivalent but more 
convenient representation of Hq reads 



H = 



tt21 



tt\2 
Hq 



where the four nodal sub-Hamiltonians Hg have been as- 
sembled into a sing le block Hi 2 , and H 21 = H^ . The 
explicit form of Hq 2 in the four-component nodal space 
is 



tt12 
Hq 



with 



g(2) 



= l~ 1 d 1 



id 2 



(4) 



where X — V or X = A. The parameter £ is a corre- 
lation length which we assume to be much larger than 
the lattice constant a. (This assumption is motivated by 
the fact that impurity scattering in real d-w&ve supercon- 
ductors seems to be predominantly forward [p3| p5|. To 
reproduce that feature, we have to take Vi mp and Aj mp to 
be slowly varying.) The strength of the disorder is mea- 
sured by the constant <?, and / is a correlation function 
normalized to unity: J cPr/(|r|/£) = 1. 

Later, we will argue that taking Vimp and Aj mp from 
the same distribution is simply a matter of technical con- 
venience. Allowing the distribution of disorder in the two 
channels to be different will not change the universality 
class, and therefore not alter the qualitative behaviour. 




FIG. 3. Impurity scattering between the nodes of 
a d-wave superconductor. The ellipsoidal lobes indi- 
cate the momentum support of the four low-energy sec- 
tors 1,1,2,2. By way of example, the four different 
scattering channels w pq with mean momentum transfer 
(p,q) — a~ [(0, 0), (7r, 0), (0, — tt), (tt, — 7r)] acting on node 1 
are indicated. 
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The stochastic part of the quasi-particle Hamiltonian 
then has the form 



^imp = Vi mp (J2 + Ai mp (7i — | JJ 12 ^ 



ttYI 
imp 



imp 



where H^ 2 p — A imp — iVi mp . Projection of this operator 
on the space of four low-energy sectors gives 



Hi 



WOO 


tU— 7T7T 


W-ttO 


WOtt \ 




tviV — TT 


WOO 


Wo-tt 


WnO 




W-xO 


WOtt 


woo 


Wirir , 


1 i 


W0-7T 


W-ttO 


W — — 


Woo 1 





(5) 



where the matrix elements are random functions of posi- 
tion given by 



u pq 



— e i(px+qy) (a. - j\A ) 

— c V imp "impj- 



They describe stochastic scattering between nodes that 
differ by a momentum (p, q) (c.f. Fig. ^). Writing 



g pq = J d r(w pq {0)w pq (r)) 

for the coupling between the nodes, and assuming the 
Gaussian disorder specified above, we have 



9pq — 9—p—q 



Note two consequences which now follow from the in- 
equality £ ^> a: (i) the variances g pq for large momentum 
transfer (p, q) are suppressed as compared to the variance 
g in the soft channel, and (ii) the various scattering chan- 
nels are approximately statistically independent. 

The sum of the pure (|]) and disordered (|J) operators 
constitutes the full Gorkov Hamiltonian: 



H = 



H 1 



H 12 = Hi 2 



H 



12 
imp 



(6) 



which forms the basis for our analysis below. Before turn- 
ing to the generic symmetries of this Hamiltonian, we 
introduce a perturbation which is both interesting and 
physically relevant. 

In the present paper we do not touch upon the sub- 
ject of perturbations involving the electron spin, but we 
will repeatedly comment on the role of time-reversal in- 
variance and the effects of breaking it. Time-reversal 
symmetry is broken in the mixed state of type-II su- 
perconductors, where a magnetic field penetrates in the 
form of vortices. Exactly how vortices influence the low- 
energy density of quasi-particle states of a pure rf-wave 
superconductor has been the subject of an ongoing de- 
bate p^-pl]|. Progress was recently made by Franz and 
Tesanovic (FT) ^8|, who suggested to perform a singu- 
lar gauge transformation similar to what is done in the 
Chern-Simons gauge theory of the fractional quantum 
Hall effect. Because each vortex carries half a magnetic 



flux quantum, a unit cell consisting of two vortices is 
used. Dividing the vortex lattice into A and B sublat- 
tices, FT make a singular gauge transformation (centered 
around vortices) in the particle sector for the A sublat- 
tice, and in the hole sector for the B sublattice. The 
purpose of the transformation is to cancel the magnetic 
field on average, thereby allowing to describe the quasi- 
particle states as Bloch waves. Although FT focus on the 
case of a perfectly regular vortex lattice, the same trans- 
formation can be applied when the vortex positions are 
given a random component. After linearization and pro- 
jection on a single node (and after rotation by e i7r<Tl / 4 ), 
the gauge-transformed Hamiltonian for, say, node 1 is: 



H 



Ho 1 



(v 



+ (T 2 {vf -«f)+7 1(J l( v l 



Here v A ' B are the supercurrent velocities for the A and 
B sublattices. We see that the differences of the veloc- 
ities add to Vimp and Ai mp , while their sum enters as a 
scalar potential (proportional to the unit matrix cr = 1% 
in particle- hole space). The latter fact was first under- 
stood by Volovik |j61|] , and is the message to carry away 
from the above discussion: the main consequence of time- 
reversal symmetry breaking by vortices is the addition of 
a scalar potential to the single-node Dirac Hamiltonian. 

Anticipating the discussion of the next subsection, we 
remark that the different treatment of the particle and 
hole sectors by FT breaks the symplectic symmetry (Q) 
of the Gorkov Hamiltonian. For our purposes, this causes 
no problems. The only agent we will need is the scalar 
potential proportional to v£ + , and this can be pro- 
duced also from a gauge transformation that treats par- 
ticles and holes on an equal footing. 



B. Generic symmetries 

An essential pre-requisite to the construction of the 
low-energy effective theory for dirty d-wave supercon- 
ductors, just as for disordered metals, is a solid under- 
standing of the fundamental symmetries of the quasi- 
particle Hamiltonian. In recent years it has become 
clear that, because of quantum interference channels 
that owe their existence to the particle-hole degree of 
freedom of the Gorkov equations, the low-energy quasi- 
particles of disordered superconductors transcend the es- 
tablished framework of the "threefold way" J6^| , namely 
the Wigner-Dyson classification scheme by unitary, or- 
thogonal, or symplectic symmetry. Dirty superconduc- 
tors generically fall into one of four non-standard sym- 
metry classes |34|] . Borrowing the notation from Cartan's 
table of symmetric spaces, these have been termed C, CI, 
D, and Dill. Specifically, spin-singlet superconductors 
(with conserved quasi-particle spin) belong to the classes 
C or CI, depending on whether time-reversal symmetry 
is broken or not. 
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While these are the generic symmetry classes of dirty 
superconductors, other classes appear, and complicate 
the symmetry pattern, if system-specific conservation 
laws of exact or approximate nature are present. As we 
have reviewed, the low-energy quasi-particles of rf-wave 
superconductors organize into four nodal sectors. In the 
physically relevant case of "soft" disorder, the scattering 
between nodes is suppressed as compared to the scat- 
tering within nodes. Of course, on very large time scales 
any small amount of scattering couples the nodes, leading 
to generic low-energy behaviour. On intermediate time 
scales, however, the nodes can be regarded as isolated. 
What symmetry class the disordered quasi-particles con- 
fined to a single node belong to, is not immediate from 
general principles. By inspection, we will find that the 
class is ^4111 if time reversal invariance is present, and 
A if not. While the latter is just the Cartan label for 
the "unitary" Wigner-Dyson class, an example for the 
former arc Dirac fermions in a random U(iV) gauge field 
(see e.g. @||). 

Having sketched the general picture, we now turn to 
the details. We begin by recalling that the Gorkov Hamil- 
tonian for a disordered spin-singlet superconductor sat- 
isfies the relation 



C : H = -a 2 H T G 2 



(7) 



which we refer to as "particle-hole" (ph) symmetry. Its 
physical origin is conservation of spin in conjunction with 
Fermi statistics, and it is violated if and only if spin ceases 
to be a constant of the motion. For systems where time- 
reversal invariance holds, which implies the absence of 
external fields and complex order parameter components 
such as id xy , the Hamiltonian matrix can be chosen to be 
real symmetric: H = H = H T . Because of the rotation 
by e I7rCTl / 4 made in going to equation (^), time-reversal 
symmetry here does not take its canonical form but is 
expressed by 



T: H = a 1 H 1 a 1 



(8) 



When both C and T are valid symmetries, the quasi- 
particles of the superconductor are said to be in class CI. 
When T is broken but C is still present, the symmetry 
class changes to C. 

The presence of both C and T constrains H to be of 
the form 



H = 



Z 



Z T = z 



(9) 



Comparing this with the previous section we notice that 
the Hamiltonian (^) is indeed off-diagonal but H 12 , un- 
like Z, is not symmetric. Looking even further back we 
see that, before linearization and projection on the four 
nodes, the off-diagonal block H} 2 = Ay — itij = Hj 2 
still obeyed the constraint of being a symmetric matrix. 
What is the origin of this apparent discrepancy? 



The answer is this. In writing cfjj 



5k 



T-k n +8k 



and mak- 



ing the assignment to nodes, we used different conven- 
tions in the "particle" (df) and "hole" {d{) sectors, so as 



to arrange for H in (|4|) to have entries only on the di- 
agonal. If we had taken the conventions to be the same, 
we would have obtained H 12 = Zq 

( 



Z imp with 



Z 











\ 



V 



(10) 



/ 



and 



^7T — 7T 


WOO 


U>0-7T 




W 


1V— 7T7T 




W 07r 


U>0-7T 


W-ttO 




woo 




WQtt 


woo 





(11) 



both of which are symmetric, in agreement with (^|). As 
follows from the expression for Zq, here one has to imag- 
ine that the pure Dirac Hamiltonian flips the particle- 
hole spinor between the particle state of one node, say 1, 
and the hole state of its conjugate, 1. This interpreta- 
tion, while perfectly valid and compatible with momen- 
tum conservation, is not the one commonly adopted in 
the literature. There, one follows the convention of in- 
terchanging the two hole states, so that the Dirac Hamil- 
tonian simply flips between the particle and hole states 
of one and the same node. To facilitate comparison with 
the lite ratur e, we tacitly made this change of basis in 
Section II A , The transformation between the two bases 



is effected by a matrix we denote by t\ : 



ri 




(12) 



By multiplying the particle-hole spinor with this permu- 
tation matrix in the hole sector (but not in the particle 
sector), the symmetry relations (0) and (||) are trans- 
formed into 



C 
T 



#=-((72 

H= (ax 



Ti)H T (a 2 <E> n) , 
Tx)H T {ax® Ti) . 



(13) 



As is easily verified, these relations are obeyed by the 
Gorkov Hamiltonian of equation (^). 

In summary, we distinguish between two choices of ba- 
sis, and hence between two ways of writing the Hamilto- 
nian. Both have their respective advantages and disad- 
vantages, and depending on the given context we will use 
one or the other, whichever is better suited. We write 

H e = UHU^ 

for the Hamiltonian in the canonical representation (|^) . 

There now exist two different scenarios. If the inter- 
node scattering is assumed to be strong, a low-energy 
quasi-particle prepared in a given initial configuration 
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will quickly attain a state where its wavefunction is uni- 
formly spread over the space of four nodes. In that case, 
the presence of the nodal exchange operator t\ in (|l^) is 
of no consequence, and we are back to the generic situa- 
tion where the symmetry class is C or CI depending on 
whether time-reversal invariance is broken or not. 

On the other hand, if the inter-node scattering in H lmv 
is negligibly weak, the four low-energy sectors decouple, 
and then the two individual operations in (^), both of 
which relate nodes to their conjugates, become ineffec- 
tive. What remains effective is an operation relating each 
isolated node to itself. By applying the symmetries C and 
T in sequence, 

H — * -(cr 2 <S> Ti)H T (a 2 <S> n) — ^ -<r 3 Ha 3 , 

we see two facts: (i) the product operation CT does not 
mix the low-energy sectors and (ii) if both C and T are 
good symmetries, the nodal sub-Hamiltonians are odd 
under conjugation by a 3 . (The latter follows more di- 
rectly from the basic equation @)). Given the presence of 
disorder, no further symmetries are expected. Thus, the 
quasi-particle Hamiltonian for an isolated node is con- 
strained only by 

H = -<T 3 Ha 3 . 

This relation is the defining equation of the "chiral" sym- 
metry class Alll. 

Any perturbation that breaks time-reversal invariance 
destroys T and hence CT. Under such conditions, the 
nodal sub-Hamiltonians do not have any symmetry other 
than Hermiticity. We say that Hamiltonians of the last 
type belong to class A, the standard Wigner-Dyson class 
with unitary symmetry. 

This concludes our analysis of symmetries. Wc have 
seen that, depending on the hard or soft nature of the 
scattering potential and the presence or absence of time- 
reversal invariance, the Gorkov Hamiltonian for dirty d- 
wave superconductors belongs to one of four symmetry 
classes: CI, C, ^4111, or A. These symmetries have phys- 
ical consequences (such as singular vertex corrections 
to the density of states, the thermal conductivity etc.), 
which arc best evaluated in a ficld-thcorctical formalism. 



C. Field-integral formulation 

Our goal is to compute the disorder average of the 
Green function, (G(E)) = ((E-H)- 1 ), and for this pur- 
pose we employ the machinery of supersymmetry \ p2\ . 
In that method, Green functions are generated from a 
Gaussian functional integral, 

Z\j,k] = J P We ^ T(£ - ff W(^+^) ; 

where Im E > is assumed and H , defined in (||) , is the 
Gorkov Hamiltonian projected on the four nodal regions 
of the ci-wave superconductor. 




are fields that have 2x2x4 components each, where 
S, S (x, k) denote complex (Grassmann) fields with 2x4 
components in the tensor product of ph and node space. 
While convergence of the integral requires the commuting 
components of the fields <j) and ip to be related by com- 
plex conjugation, the anticommuting components of the 
two fields are independent. As usual, Green function ma- 
trix elements are obtained by differentiating twice with 
respect to the sources j and k. 

A noteworthy feature due to the symmetry C in ( |l3| ) 
is that retarded (Im E > 0) and advanced (Im E < 0) 
Green functions are related by 

G{E) = -(<r 2 ® T 1 )G T {-E)(a 2 ® n) . 

This has the consequence that the functional Z[j, k] gen- 
erating a single Green function can also be used to com- 
pute the disorder average of the two-particle Green func- 
tion (G(0+)G(0-)) at zero energy (E = ± = ± iS). 
As the emphasis here will be on the basic structure of 
the theory (rather than on specific observables), we tem- 
porarily suppress the source content and focus on the 
functional Z[Q], 

In the next step, we are going to adapt the func- 
tional integral to the particular symmetries of the Gorkov 
Hamiltonian. In normal conductors, time-reversal sym- 
metry is known to give rise, via the so-called Cooperon 
mode, to quantum interference corrections to diffusion 
that are infrared singular in dimension d < 2. Simi- 
lar modes appear in the present case |54j, as a result of 
the discrete symmetries C and T. In the impurity dia- 
gram technique, these modes emerge as divergent series 
of ladder graphs, or maximally crossed diagrams. For 
a semiclassical interpretation of these modes we refer to 
J34| . The field-theoretic approach elevates their status to 
those of Goldstone modes due to the breaking of certain 
global continous symmetries, which in turn derive from 
C and T. 

In view of this, we now exercise special care to trans- 
late the fundamental symmetries of the Hamiltonian into 
symmetries of the functional integral. For that purpose 
we find it convenient to transform to the representation 
UHU^ = H c = ciRcZ — <Ji\\aZ introduced in the previ- 
ous subsection. Using it, we write 

Z[0] = J cxp (i£(/3 T 7 + /3 T 7) - */3 T Z 7 - iP T h) , 

where the unitary transformation U was absorbed into 
the measure of the functional integral (which is implicit, 
as is integration of the exponent over position space). 
/3,/3 (7,7) are the PH components of (j) {ip). 
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Now recall that the operators Z, Z are symmetric. In 
order to firmly install this property in the functional in- 
tegral, we symmetrize and "double" ]52| the field space: 

/3 T Z 7 = i/3 T Z 7 + ZP ee * S Z* . 

Here the matrix a^ F acts in the boson-fermion space of 
commuting and anticommuting field components, and 
the spinor and its transpose \l/ s are defined by 



J_/ 7 



1 



(14) 



We refer to the two-dimensional space comprising the 
upper (7) and lower (/?) components of the spinor \P as 
"charge conjugation" (cc) space. Note that the scalar 
product determined by the transposition rule W 1— > V E' S is 
symmetric: \I/ S SI/' = ty'ty. Doing the same symmetriza- 
tion for the other terms in the exponent, we obtain 



Z[0]= / V^V^c 



- fed 2 



£ = £ 



D 



(15) 



A node-resolved representation of the Lagrangian is got- 
ten by making the decompositions 



* EE 




# = 




The overbar here does not mean complex conjugation. 
Following the conventions of conformal field theory, we 
use it to denote fields whose correlation functions are an- 
tiholomorphic in the limit of zero disorder. In the final 



step, we substitute the nodal decomposition into (15) 
which gives 



r 12 



r 1 



c 2 



f 12 
>-T> 



2^{(vd {1) + iw 00 )^i + 2^ s l (vd (1) + iw 00 )^i 
+iiu 7r _ 7r *J* 1 + f * J 

+iw n -^i$x + iw-^{ *i , 

2%(vd (2) + iw 00 )^2 + 2%{vd {2) + mD o)*2 
+iw_ 7r _ x *2*2 + M«W*f*2 
+iw)_ T _ w *2*2 + i«W*§*2 ' 

: 2^0-^*1*2+^0*1*2 
+ W_ 7r0 *f*2 + W 07 r*f*2 
+ + ^0*1*2 

+ ^0*1*2 + W 07r *?*2) > ( 16 ) 



where and C 1 ^ govern the pairs of nodes (1,1), 

(2,2), and the coupling between them, respectively. The 
factors of two in these expressions arise from combining 
terms: V s A$ ± ^> S A$>' = 2^' S A^ for a symmetric (or 
antisymmetric) operator A. 

Equations (JlJ) and ( |l6| ) cast the dirty d-wave problem 
into the field-theoretical form that all subsequent analy- 
sis will be based on. 



D. Symmetries of the Gaussian field theory 



We next explore the symmetries of the Lagrangian 
(|l6|). The outcome will pre-determine the structure of 
the low-energy effective field theories that derive from it. 

To prepare the stage, we recall that the word "symme- 
tries" has two aspects to it. Firstly, there exists a symme- 
try group, which acts by transformations that leave the 
Lagrangian invariant. This is the proper meaning of the 
word "symmetry" . Secondly, the degrees of freedom of 
the theory, the fields, take values in a target manifold. In 
current speak the latter, too, is sometimes referred to as 
the "symmmetry" of the theory, although it is of course 
different from the symmetry group. The remark we wish 
to make is the trivial and yet important statement that 
these two meanings must be kept apart. In the present 
section we are going to discuss the first aspect, namely 
the symmetry group of the action functional. The struc- 
ture of the target manifold will be discussed in Section 
[II, after the introduction of the relevant low-energy de- 
grees of freedom. 

We begin by considering the case of highest symme- 
try where all of the disorder except for woo is switched 
off. Symmetry breaking caused by hard scattering w pq 
and the perturbations that break time-reversal invari- 
ance, will be discussed afterwards. 

Class yffll (forward scattering, T-invariance) : on set- 
ting E — and w pq = (for pq ^ 00), the Lagrangian 
reduces to four independent terms: 

L = 2^ 6 1 {vd ( - 1) + iw o)*i + 2*? (vd w + iw 00 )*i 
+ 2*2>d (2) + zw 00 )*2 + 2§ 2 : 0<9 (2) + ™ 00 )*2 , 



which are distinguished by the dichotomy of holomor- 
phic (*i>) versus antiholomorphic (*£) fields, and by the 
grouping into pairs of nodes: (1, 1) and (2,2). This La- 
grangian describes (anisotropic) Dirac fermions in a ran- 
dom Abelian vector potential A — wqo, A — wqo. Since 
the four terms are identical in structure, but involve dif- 
ferent partial derivatives eK 1 -* 7^ for non-vanishing 
anisotropy 7^1, the symmetry group will be a Cartesian 
product of four copies of the same group, which is read- 
ily identified as the supergroup GL(2|2). Indeed, looking 
at, say, the first two terms, and noting that u>oo(r) and 
wqq(t) are just complex numbers, we see that a local 
transformation 

¥f(r) 1 * ^(r)r(z!)- 1 , *!(r) 1 * T(*i)*i(r) , 
*f(r) ~ ^(r)T(f!)- 1 , %(r)^T(2i)*i(r) , (17) 



where z\ = 7x1 + 1x2 and z\ = 7x1 — 1x2, leaves the ac- 
tion functional invariant. Because ^1,^1 comprise two 
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commuting and two anticommuting components, the in- 
vertible supermatrix T is of size (2 + 2) x (2 + 2). Hence 
the symmetry group of the present case is GL(2|2) or, 
rather, four independent copies thereof. As we shall see, 
the local nature of the symmetry group makes the present 
theory exactly solvable. Note that here, and through- 
out this section, the word "symmetry group" means the 
complexified |6q | symmetry group, which ignores the re- 
lation that exists between the bosonic field components 
and their complex conjugates. 

We mention in passing that the symmetry group of 
the totally clean limit (u>oo = woo = 0) described by free 
fields, is four copies of the orthosymplectic supergroup 
OSp(4|4). This will come to play a role when we turn to 
the method of non-Abelian bosonization. 

Class A (s oft scattering, no T-invariance) : as we saw in 
Section II A , the Hamiltonian H e ceases to be off-diagonal 
in ph space (while remaining diagonal in CC and node 
space) when time-reversal symmetry is broken. Hence 
the breaking of T introduces terms into the Lagrangian 
that mix the two sectors of holomorphic and antiholomor- 
phic fields. Consequently, invariancc of the Lagrangian 
requires the actions of GL(2|2) in these sectors to be re- 
lated to each other. This "locking" of group actions is 
the only effect T-breaking has on the symmetries, and 
therefore the net result is that the four copies of GL(2|2) 
get reduced to two copies, corresponding to the two pairs 
of nodes (1,1) and (2,2). 

Let us look at the locking of GL(2|2) grou p actions in 
more detail. As we recalled in Section II A , the super- 
current flow in the mixed state acts as a scalar potential 
(v A + v B )l ® r 3 . After transformation to the canonical 
representation, H e = UHU\ this reads (v A + i> b )0 3 ®t 3 . 
Consider therefore the T-breaking perturbation 

ip T (o- 3 <g> r 3 )0 = /3 T r 37 - /3 t t 3 7 = 2* s (cr 3 3C ® r 3 )* . 

In order for this perturbation to remain unchanged un- 
der the transformation (|l7|), we must require T s = rjT -1 !] 
where 77 = cr 3 c ® r 3 . We will use this relation to identify 
the target manifold for class A in Section [II. 

Class CI (hard scattering, T-invariance) : we next con- 
sider the situation where the hard scattering channels 
are switched on. The matrix elements w pq ^ then 
mix fields pertaining to different nodes, whence the four 
copies of GL(2|2) for class AIII no longer act indepen- 
dently but become locked to each other. Moreover, in 
order for terms in the Lagrangian such as 5 , f^[ , 1 to be 
invariant under $i(r) t— > T<I>i(r), we require T S T = 
1, where the supermatrix transpose T s is defined by 
(T$) s = $ S T S . From Q we read off that $ s is ob- 
tained from <I> by taking the ordinary transpose, and then 
exchanging the two commuting (anticommuting) com- 
ponents by multiplication with o~\ (io~2)- The equation 
g T <J\g — <j\ defines an orthogonal group 0(2, C), and 
5 T o"25 = 02 the symplectic group Sp(2,C). Thus the 
condition T S T = 1 determines an orthosymplectic sub- 
group OSp(2|2) of GL(2|2). The full symmetry group 



consists of two copies of OSp(2|2), since the holomor- 
phic fields of the pair (1, T) and antiholomorphic fields of 
(2, 2) can be transformed independently of the remaining 
fields. 

The following remark may be helpful. Given the phys- 
ical distinction between the BB and ff sectors, there are 
two versions of OSp to consider. The first one acts by 
symplectic transformations in the former sector and by 
orthogonal transformations in the latter. (We refer to 
this as "symplectic" bosons and "orthogonal" fermions 
for short.) The OSp(4|4) symmetry group of the free- 
field problem is of this kind. In the second version of 
OSp, which is the relevant one here, the roles of the or- 
thogonal and symplectic groups are interchanged ("or- 
thogonal" bosons and "symplectic" fermions). 

As another aside, we mention that en route from soft 
to hard scattering one could imagine a scenario where 
the scattering couples each node only with its conjugate 
§. Such a scenario would lead to a different symmetry 
group and eventually to a different target manifold [jl9| . 
However, we see no physical room for this theoretical 
possibility and will not pursue it here. 

Class C (hard scattering, no T-invariance): the sym- 
metry group for this final case can be approached from 
class A by including hard scattering, or from class CI by 
breaking T. In either way, one finds that the remain- 
ing symmetry is just a single copy of the group OSp(2|2) 
with orthogonal bosons and symplectic fermions. This is 
the orth osymp lectic supersymmetry omnipresent in class 

For future reference, the symmetry groups for the four 
different cases are summarized in the following table, 
where the notation x"G means n independent copies of 
the group G. 



T 


hard scattering 


symmetry group 


+ 
+ 


+ 
+ 


x 4 GL(2|2) 
x 2 OS P (2|2) 
x 2 GL(2|2) 
OS P (2|2) 



E. Disorder average 

We now carry out the disorder average. By integrating 
the Gaussian generating functional over the distribution 
of the matrix elements w pq , which are subject to the cor- 
relation laws stated in Section [II A| , we obtain 

Z av [0] = (Z[0]) = J V^V^e-I {CD+CE)d2r , 

where 

£ D = 2vC + £ dis , (18) 
and all perturbations O pq are quartic in the fields: 
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Goo = 4*;*!*;*! 



44-^2*2*2 > 
4\l/j\l/ I \['2*2 ' 



-Sab 



+ 4*f* 2 *?*2 + 4*1*2*1*2 > 
+ *^*2*i*2 + *2*2*2*2 • 

In the initial stage of the calculation, disorder averaging 
produces expressions of the form 



d 2 r 



dV/(|r-r'|/0 (* s *)(r) (* s *)(r') 



where / is the correlation function of the disorder. To 
arrive at the above expressions for O pq , we omitted the 
finite spread of this function due to a non-zero value of £. 
This approximation is justified by the fact that all cor- 
rections from expansion around the local limit r — r' = 
carry at least two derivatives, which renders them irrel- 
evant in the renormalization group (RG) sense. Note 
that all of the perturbations O pq are marginal by power 
counting at the free-fermion point. 

Eq. ([L8ft defines the disorder-averaged theory. Before 
going further, we will inquire into the nature of the RG 
flow caused by the perturbations O pq . 



F. Renormalization group 

In this section we perform a one-loop renormalization 
group analysis to explore the relevance of the five cou- 
plings g pq . There exists a standard formula |6^|3!| to 
use for that purpose, which refers to the operator prod- 
uct expansion (OPE) of the perturbating operators. The 
general statement is that, if the short-distance expansion 
for a set of marginal perturbations O^' has the form 



O {t \r)O^(0)d 2 r = 2njJ2ckO {k) (V) 

k 



t<\r\<l+St 



the corresponding couplings gu renormalize according to 
the equation 

Me) = -nj^c'k9i9j , 



d\n£ 



ij 



where I is the short-distance cutoff, /?/. are the beta func- 
tions, and c l j! are called the structure constants of the 
algebra of operators O' 1 '. 

To apply this formula to the problem at hand, we need 
the OPE for the fundamental fields. Apart from an over- 
all multiplicative constant, which can be removed by a 
conformal rescaling of the fields, these are 



#1^(21, ^2)*i,b(0,0) = 



Sab 



aX2 


+ ia 


1 x 1 




Sab 




ax i 


— ia" 






-Sab 




axi 


+ ia" 


X X 2 



*M(a:i,£2)*2,.B(0,0) = 



Here a = = \/t/ A is the (square root of the) 
anisotropy parameter, and A, B is a composite index 
built from the BF and CC indices of the fields. As usual in 
this context, the above relations have to be understood 
as identities that hold under the functional integral sign. 

By using the free-field expansions in conjunction with 
Wick's theorem, we can now work out the OPE for the 
set of composite operators O pq . By a straightforward if 
tedious calculation, this yields the RG flow equations 



9 
9 

?7T0 



0, 

2^7T7r 
33tt0 

2(55*- 



(19) 



ax2 — ia L X\ 



where gi stands for dgi/dhii (times some unimportant 
constant which we do not specify). Notice the following 
features: 

• The beta functions contain no terms of linear order 
in the couplings, which expresses the fact that we 
are dealing with a set of marginal perturbations. 

• In the physical regime of positive couplings, the 
beta functions are never negative, so none of the 
couplings decreases under the RG flow. 

• None of the couplings supports its own flow. In 
particular, for the model with only forward scat- 
tering (g pq = 0) the coupling constant g is truly 
marginal. 

• For generic initial data, the nature of the flow is 
marginally relevant, i.e. the couplings (save for the 
constant g') increase under the flow, although the 
rate of increase vanishes in the limit of weak disor- 
der gi \ 0. 

• The anisotropy parameter 7 does not enter the one- 
loop RG equations. 

The fact that the anisotropy has disappeared from the 
RG equations can be understood heuristically as follows. 
The OPE scheme for the set {O pq } encodes the one-loop 
renormalization of the four-fermion vertices of the the- 
ory. These vertices are related to the fermion self energy 
through a Ward identity. To explore the effects of the 
anisotropy, one may therefore directly analyse the self- 
energy diagrams. The latter have been shown || to be 
unaffected by the anisotropy, to leading (or one-loop) or- 
der. The physical reason is that the first-order self en- 
ergy diagram measures the Born scattering rate between 
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the four low-energy sectors. This rate is not affected by 
the ellipsoidal shape (expressing the anisotropy) of the 
low-energy lobes in momentum space; it only depends 
on the total phase volume available. In higher orders of 
perturbation theory, the situation changes and the self 
energy begins to be affected by the anisotropy in mo- 
mentum space. In particular, higher-order scattering be- 
tween neighbouring nodes is suppressed || . Accordingly, 
we expect that the beta functions at higher loop orders 
do depend on the anisotropy of the model. 

Let us stress two messages that emerge from the cur- 
rent section: (i) owing to the quadratic dependence of 
the beta functions on their arguments, the RG flow is 
marginal (or very slow) in the limit of small couplings 
(or weak disorder), and (ii) for generic values of the 
disorder-generated couplings, the renormalization group 
flow drives the system away from the free-fermion the- 
ory. In combination with the phenomenological input 
that the inter-node couplings g pq are small as compared 
to the intra-node coupling g, the first message leads us 
to expect a prolonged crossover intervening between a 
ballistic regime at short scales, and the diffusive regime 
governed by a non-linear sigma model at large scales. 



III. HOW TO PROCEED? 

Anticipating that a direct perturbative analysis of the 
theory ( |l8f ) will not capture all of the important physics, 
one is tempted to employ the "standard" scheme for 
dealing with weakly disordered non-interacting parti- 
cles: Hubbard-Stratonovich transformation followed by 
a saddle-point analysis and gradient expansion. It turns 
out, however, that this approach is beset with a num- 
ber of problems, none of which appears in systems with 
a non-relativistic kinetic energy. These difficulties will 
eventually force us to adopt a strategy where elements 
of the standard scheme are supplemented, and even su- 
perseded, by an alternative scheme tailored to relativistic 
fcrmions: non-Abelian bosonization. 

For pedagogical reasons, we will here proceed in a con- 
servative fashion and continue somewhat further along 
the much-trodden path (Hubbard-Stratonovich transfor- 
mation and so on), introducing along the way a few con- 
cepts of general validity. In particular, we will identify 
the degrees of freedom of the hierarchy of low-energy ef- 
fective theories, and the way they are asso ciate d with 
the global symmetries discussed in Section II D| . As a 
by-product, we will be able to make the connection to 
previous self-consistent approaches. 

To begin with, we set the couplings g pq to zero and 
focus on the functional integral 



Z ft [0] = Z av [0] 



'Jvq 



containing only the largest perturbation, gOoo + s'Ooo- 
The first step of the standard approach is to make a Gaus- 
sian transformation, called the "Hubbard-Stratonovich 



transformation" . Introducing two auxiliary fields P and 
Q, and exploiting the equality of the unrenormalized cou- 
pling constants g — g' , it is straightforward to show that 

Z is [0] = J V^V^ J VQVPe'vf STr ( ( 3 2 + p2 ) 

xe - J (£ E +2 U £ +* s (iQ-P)*+* s r 1 (iQ+P)T 1 *) 



x exp -iSTrln 



E 

i 



iQ — P iZo 
iZ Tl {iQ + P) Tl 



where the values of P and Q are supermatrices of linear 
size 4x2x2 acting in nodal space and internal space (the 
latter being the tensor product of charge conjugation and 
boson- fermion space). For simplicity we here ignore the 
issue of convergence of integrals, deferring a more strict 
treatment until Section IV. 



We mention in passing that it can be understood at this 
stage why the restriction to equally distributed disorder 
in the normal and order parameter channels makes no 
essential difference: one can relax this assumption by in- 
troducing two more Hubbard-Stratonovich fields, R and 
S, coupling to the PH matrices u\ and 02-, respectively. 
On varying the action with respect to these fields, one 
finds that they vanish on the saddle-point level. Since 
a non- vanishing saddle point is a necessary condition for 
the formation of Goldstone modes, which control the in- 
frared behaviour, we conclude that these fields and hence 
the choice of (un)equally distributed disorder, are of no 
importance. 

With the aim of subjecting the functional integral to 
a stationary-phase analysis, we now vary the action with 
respect to the fields Q and P to generate a set of saddle- 
point equations. We first look for solutions (Pq, Qo) that 
are diagonal in both nodal and internal space. For van- 
ishing energy E, the solutions are 

Po = 0, 
Qo = -i« , 

where the real parameter k is defined implicitly through 

2 _ f d 2 k 1 
~g~ J (2^)2 /c 2 + e(fc) 2 ' 



(20) 



This equation, which can be regarded as the Dirac analog 
of the SCBA equation for disordered normal-conducting 
systems (with —Qo representing the self energy), has 
been discussed in the literature (see, e.g., Refs. |@,||). 
Its solution is 

k 2 = M 2 e- 8 " At /s , 

where fj, represents a UV cutoff that regularizes the inte- 
gral over momenta, which would otherwise diverge loga- 
rithmically. The parameter k can be regarded as a self- 
consistently determined self energy. This interpretation 
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is strengthened by re- instating sources j, k into Z& [0] , 
and then using this generating functional to compute ma- 
trix elements of the physical Green function. Doing so, 
one finds that these are given by the inverse of the oper- 
ator 



k — iE 
iZ 



iZ 
k — iE 



(21) 



and the information content of the above saddle-point 
approximation coincides with that of the diagrammatic 
SCBA approach. 

However, as we know from experience with disordered 
normal-conducting systems, there exists much physics 
that does not unfold at the saddle-point level, but re- 
sides in the fluctuations of Q. The same is true here. 
The saddle point (Qo,Pq) — (— IK, 0) breaks the symme- 
try x 4 GL(2|2) for class yllll as recorded in the table at 
the end of section [ID. By the Goldstone mechanism, 



this leads to the appearance of massless modes, and one 
expects that it is the fluctuations in these "soft" modes, 
rather than any single saddle point, that determines the 
behaviour of the system at large scales or low energies. 

Let us therefore ask what happens when the functional 
integral Zf 8 [0] is subjected to a global symmetry transfor- 
mation taken from x 4 GL(2|2). (As we saw, the symme- 
try is actually local, but we should hesitate to draw any 
conclusions from that, as the axial part of the symmetry 
is anomalous; see Appendix |A 1| ). Since these transfor- 
mations do not mix the nodal sectors (1, 1) and (2, 2), we 
can concentrate on one of them, say (1, 1). Consider the 
action functional evaluated at the saddle point, 



-STrln ( . 

2 I ivT 2 d (l > 



-IVT2 
K 



dC 1 



where t 2 denotes the second Pauli matrix acting in nodal 
space, and the energy E has temporarily been set to zero. 
Then recall the transformation with the matrices 
(T, T) £ GL(2|2) x GL(2|2), which are now taken to be 
constant in space. By transferring this transformation to 
the argument of the logarithm, we obtain the "rotated" 
action 
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where M — TT S . This expression states that not only 
the diagonal matrix in, but in fact any GL(2|2) config- 
uration in x diag(M s , M _1 ,M, (Af s ) _1 ) is a solution of 
the saddle-point equation. Put differently, the saddle- 
point manifold for the nodal sector (1,1) is isomorphic 
to GL(2|2). On including an identical factor for the other 
sector, (2, 2), the total saddle-point manifold of the model 
with only forward scattering becomes GL(2|2) x GL(2|2). 
(Note that a complete description of the saddle-point ap- 
proximation would have to specify the real submanifold 



to be integrated over. The present discussion gives only 
the complex saddle-point manifold for simplicity.) 

We note in passing that what we are encountering here 
is reminiscent of the phenomenon of chiral symmetry 
breaking in quantum chromodynamics. The appearance 
of a chiral symmetry group Gl x Gr = GL(2|2) x GL(2|2) 
(still maintaining the focus on a single nodal sector), with 
independent left and right factors, is directly connected 
with the oddness of the Hamiltonian under conjugation 
by (T3 in ph space: o~%Ho~?, = —H. The latter symme- 
try is broken by the saddle point Qq = —iK, which is 
even under conjugation by 03. In four space-time dimen- 
sions, the role of 1173 = o~\g<x is played by 175 ~ 70717273, 
for massless fermions one has 75^/75 = — H, and M be- 
comes the pion field. In QCD, as in the present case, the 
saddle-point manifold is given by a diagonal subgroup 
Gl = Gr (the "axial" symmetry transformations) of the 
chiral symmetry group. The main difference turns out to 
be that chiral symmetry is not truly broken (although its 
axial part is anomalous) in the vacuum of our 2d theory. 

For the case of forward scattering with T-invariance 
(class AllI) the above discussion a nswe rs the question, 
raised at the beginning of Section II D, concerning the 
relationship between the global symmetries of the field 
theory and its degrees of freedom. The former are repre- 
sented by T and T, and the latter by M = TT S . The next 
step of the standard scheme would be to allow M to vary 
slowly and carry out a gradient expansion of the action 
functional. However, for reasons that are spelled out be- 
low, we do not pursue this approach here but will switch 
to an altogether different strategy in the next section. 
Before doing so, we extend the construction of saddle- 
point manifolds to the other symmetry classes: CI, A, 
and C. The procedure is always the same: we let the 
global symmetry group act on the diagonal saddle-point 
Qq = —in, and the saddle-point manifold is then simply 
the result of this group action. 

Class CI (hard scattering, T-invariance): Recall from 
Section |ll D| that the couplings g pq , which mix the nodes 
and lead from class AllI to CI, reduce the symme- 
try group from four copies of GL(2|2) to two copies of 
OSp(2|2). Because there are no other changes, the above 
form of the symmetry-transformed action functional con- 
tinues to hold, with M still given by M = TT S _. The or- 
thosymplectic group property, T s — T^ 1 and T s = T , 
entails M 8 = M~ x . As T and f vary over OSp(2|2), so 
does M, and hence the (complex) saddle-point manifold 
is isomorphic to that group. The symmetry group is still 
chiral, acting on the field M(r) independently on the left 
and right by M(r) i-> TM(r)f s . 

Class A (soft scattering, no T-invariance): Again, we 
focus without loss on the pair (1,1). According to the 
discussion in Section II D . breaking of T-invariance by 



a supercurrent flow locks the left and right actions of 
GL(2|2) to each other by f s = rjT^-q with 77 = of' ®t 3 . 
Because M and T are node-diagonal, this equation sim- 
plifies to f s = ofT^of. The expression for M = TT S 
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T/i by elements h that leave Cg 
These are easily seen to form 
a subgroup GL(1|1) x GL(1|1) of GL(2|2). Hence the 
complex saddle-point manifold is isomorphic to 

GL(2|2)/(GL(1|1) x GL(1|1)) . 

This coset space is very familiar plq ] as the target man- 
ifold of the non-linear sigma model for systems in the 
Wigner-Dyson class with unitary symmetry (alias class 
A). The standard symbol for the field in this case is 
Q = Ma^ c . Identical statements apply to the pair (2, 2). 

Class C (hard scattering, no T-invariance): As before, 
the constraint T s = a^ c T~ 1 a^ c from T-breaking leads to 
Q = TafT- 1 . What is different now is that the field 
is completely locked in node space, and T £ OSp(2|2) 
obeys the condition T s — T _1 . Again, the field Q does 
not change when T is multiplied on the right by an ele- 
ment h that fixes a^ c under conjugation. Such elements 
turn out to form a subgroup GL(1|1), and Q therefore 
parametrizes the coset space OSp(2|2)/GL(l|l). 

This concludes our discussion of saddle-point mani- 
folds. A summary of the essential results was tabulated 
in the introductory section. Let us mention that major 
elements of that table have appeared in the literature, 
although in the framework of the fermion-replica trick: 

• Prior to the celebrated discovery that the order pa- 
rameter of high-temperature superconductors has 
d x 2_ y 2 symmetry, Oppermann |35| had studied a 
system he called the "superconducting Ising glass" , 
and derived for it a non-linear sigma model over 
the compact group Sp(2r), with r = owing to 
the use of the replica trick. Oppermann's super- 
conducting Ising glass is a time-reversal invariant 
system with local superconducting order and con- 
served spin. Thus it has all the prerogatives of class 
CI, and the target manifold Sp(2r) is just what one 
expects on the basis of the general classification 
[ p0| . For the important case of <i-wave supercon- 
ductivity, the Sp(2r) non-linear sigma model was 
recently rediscovered by Senthil et al. |g] . 

• According to Nersesyan, Tsvelik, and Wenger ||, 
the effect of impurity scattering on an isolated node 
of a d-w&ve superconductor is described by a non- 
linear sigma model (more precisely, by a WZW 
model) on U(r). We will discuss a supersymmetric 
variant of that model in some detail later. 

• The non-linear sigma model with target manifold 
Sp(2r)/U(r) was identified as the low-energy ef- 
fective theory for dirty <i-wave superconductors in 
class C by Senthil et al. Q. The supersymmet- 
ric extension of that theory had appeared in ]7C[] 
(to describe the physics of low-energy quasi-particle 
excitations in disordered superconductor/normal 
metal junctions) and |56| (as a description of quasi- 
particles in the core of a disordered vortex.) 



The appearance of the Wigner-Dyson class of uni- 
tary symmetry (class A) for c£-wave superconduc- 
tors in the mixed state is implicit to several recent 
papers |lp|||. 



The merit of the present work is that it assembles the 
various symmetry classes and field theories into a single 
coherent scheme. 

After this extensive tour of symmetries and their field- 
theoretical realization, we are now ready to attack the 
concrete goal of giving a bona fide construction of the 
low-energy effective theories. Equipped with a strong 
background in the supersymmetric theory of disordered 
metals, one would deem the strategy to follow quite ob- 
vious: one should decompose the Hubbard-Stratonovich 
fields P, Q into Goldstone (or massless) modes M and a 
complementary set of massive modes. One would then 
integrate over the latter in Gaussian approximation and, 
finally, expand to lowest order in the gradients VM, mea- 
suring the energy cost due to spatial variations in the 
massless fields. It came as a surprise to us that this pro- 
gram is corrupted by two sources of severe difficulty: 

• In the course of carrying out the gradient expansion 
we are confronted with a two-dimensional version 



of the chiral anomaly (see Appendix A 1 ) 



• For the physically interesting limit of systems with 
T-invariance and forward scattering only (class 
AIII), we are unable to isolate and eliminate the 
massive fields in a controlled manner. 

The presence of an anomaly leads to disastrous results 
when the gradient expansion is carried out in the most 
straightforward way. This, however, is a surmountable 
problem; once the nature of the anomaly has been un- 
derstood, a less naive and properly regularized expansion 
scheme can be set up to yield the correct answer. The 
second problem is more serious. Its origin is that Dirac 
fermions are perturbed by a random Abelian vector po- 
tential (which is the concrete form taken by the disor- 
der in the present realization of class Alll) in a truly 
marginal way. What this means is that the system re- 
tains the conformal invariance of the free-fermion theory, 
with no mass scale being generated. As a result, any 
attempt to integrate out the "massive" modes perturba- 
tively is doomed to fail. In principle, one could try to go 
beyond the Gaussian saddle-point approximation for the 
massive modes but, in practice, this seems unfeasible. 

In this situation it comes as a relief that an alterna- 
tive approach to the problem exists: the method of non- 
Abelian bosonization due to Witten [B8|. 



IV. NON- ABELIAN BOSONIZATION 

The utility of non- Abelian bosonization as a tool to 
construct the low-energy effective action for the quasi- 
particles of a d-w&ve superconductor was first recognized 
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by Nersesyan, Tsvelik, and Wenger ||. These authors 
used the fermion-replica trick to deal with Dirac fermions 
perturbed by various types of disorder. By bosoniz- 
ing the disorder-averaged fermion-replica theory, they ar- 
rived at a Wess-Zumino-Novikov-Witten (WZW) model. 
This model had the attractive feature of exact solvability, 
which enabled NTW to predict some exact values for the 
scaling exponent of the low-energy density of states. As 
discussed in the introduction, however, the original ap- 
proach of NTW has a number of weak spots, one of which, 
namely the neglect of the fact that the multi- valued term 
in the WZW action changes sign when the orientation of 
position space is reversed, bears drastic consequences. 
Following NTWs original strategy, it is the goal of this 
section to introduce a bosonization scheme in which (i) all 
symmetries of the system are fully included, (ii) the de- 
pendence of the WZW action on orientation is accounted 
for, and (iii) the physics of the crossover regimes in be- 
tween the limiting cases of hard and soft scattering is in- 
cluded. We will begin by reviewing a few elements of the 
general bosonization approach to supersymmetric theo- 
ries, and then discuss the application of these concepts 
to the case at hand. 

Originally, non-Abelian bosonization was introduced 
to construct bosonic representations of relativistic two- 
dimensional fermion models with continuous internal 
symmetrie s |p8| . More recently, the approach has been 
extended [|49f to the supersymmetric case of relativistic 
fermions supplemented by a bosonic ghost system, which 
is the setup needed to treat problems with disorder. Be- 
low we will review some elements of this approach, in 
a form tailored to the d-wave application. Readers who 
are familiar with bosonization will find that there are no 
structural differences to the original, fermionic version. 
Readers who are not may wish to consult first an intro- 
duction to standard non-Abelian bosonization, e.g. Wit- 
ten's original and highly pedagogical article |3q| . A de- 
tailed account of the principal supersymmetric extension 
of the approach can be found in Ref. (49| . 

Before turning to the problem of the d-wave supercon- 
ductor, let us begin with some preliminary considerations 
of a generic model: consider the functional integral 

Z = J P^e- S °^](...)^-, 

where ip and tp are fields with 2m bosonic and 2m 
fermionic components and the ellipses (. . j, stand for 
certain operators constructed from these fields. The ac- 
tion functional So denotes the free supersymmetric action 

S [tp, ft = J d 2 r (ft dip + ft dip) , (23) 

where tp i— ► ft is an "orthosymplectic" transpose which 
differs from the s-operation introduced earlier in that it 
defines a skew symmetric scalar product: 

ft ft = -ftt^p . 



(Although skew symmetry is the only property that mat- 
ters for our purposes, it may be helpful to recall the con- 
crete realization given in |^] for n = 1: ft = (/3, 7, 6, — c) 
and ip = (7, /3, c, b) T , where Greek and Roman letters de- 
note fermionic and bosonic components, respectively.) It 
is implicitly assumed that the existence of the functional 
integral is ensured by the presence of some convergence 
generating term in the action. (For an explicit example 
of such a term, see below.) 

Furthermore, let M be a field taking values in the 
(anomalous) symmetry group of So, which is the complex 
supergroup OSp(2m|2m) of matrices fulfilling M*M = 1, 
with the matrix operation M 1— > M* defined by (Mip)* = 
ftM*. The basic statement of non-Abelian bosonization 
then is that the functional integral Z defined above is 
equivalent to 

Z'= f OTe- ff M(...) M , 

JC\D 

where 

W[M] = -JL J d 2 r STr (d.M^d.M) + ^T[M] , 

T[M] is the WZW functional, and the ellipses stand for 
a representation of (. . .),/, ^ in terms of the matrix field 
M to be specified momentarily. The phrase "equiva- 
lent" here means that, with the transcription (. . j <-> 
(. . .)m understood, the two functional integrals Z and Z' 
produce identical results. 

Following Witten, an explicit expression for T can be 
written down by choosing some extension of the field M 
to a one-parameter family of fields M(t) with the prop- 
erty M(Q) = 1 and M(l) = M. Then T[M] = 

J d 2 r dt e x ^ v STr (M^ 1 d\M M~ l d M^ 1 d v M) , 

where (A, n, v) represents the coordinate triple (t, x\, £2), 
and e A/ily is the fully antisymmetric tensor. 

The bosonization dictionary needed for the transcrip- 
tion (. . <-> (. . .)m is given as follows: 

tpft <-> J = (2w)- 1 MdM- 1 , 
tpft <-> J ee (2tt)- 1 M~ 1 BM , 
iPft <-> l~ x M , 

4>ft ^ r l M- 1 , (24) 

where I is some length scale serving to UV-regularize the 
Dirac theory. 

Finally, the subscript "C\D" in the definition of Z 1 in- 
dicates that the functional integration over M does not 
really extend over OSp(2m|2m). The reason is that the 
invariant metric STr(M~ 1 dM) 2 on this group, as on any 
complex group, is not Ricmannian. This in turn im- 
plies that a functional integral controlled by the metric 
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term STi (d fJi M~ 1 d fl M) cannot be defined by an unre- 
stricted integration over OSp(2m|2m). (This term has 
indefinite sign, which spoils the convergence of the in- 
tegral.) For ordinary groups the obvious remedy is to 
restrict the complex group to a compact real subgroup; 
one passes, for example, from 0(2m,C) to 0(2m). No 
such remedy exists for supergroups such as OSp(2m|2m), 
as these do not possess subgroups on which the invari- 
ant metric becomes of definite sign. To rescue the func- 
tional integral, one has to abandon the group structure 
and restrict OSp(2to|2to) to what is called a Rieman- 
nian symmetric superspace of type C\D in the terminol- 
ogy of Rcf. p0[ . Referring for a detailed discussion to 
Rcf. Q|, we here merely mention that the bb sector of 
this manifold is given by the non-compact coset space 
Sp(2m, C)/Sp(2m), while the ff sector is the compact 
group 0(2m). 

We next apply this general apparatus to the d-wave 
system. Firstly, it is necessary to cast the free action 
J cPrCg (we here set v = 1) of the system into the form 
of the action Sq above. To this end we define 
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try), but tpi^ = +'02'0i- O ne then verifies that the free 
part of the Lagrangian can be represented as 

n=l,2 

C E = -2iE Y ^V„ ■ 

71=1,2 

To make this Lagrangian amenable to a direct applica- 
tion of the bosonization rules, we temporarily remove the 
anisotropy inherent to the operators d^ n >. This can be 
done by rescaling the coordinates for the sector of nodes 
(1, 1) according to x\ — > r )~ 1 x\,X2 — > £2 and the coor- 
dinates for the complementary sector (2, 2) according to 
x\ — ► X\,xy, — ► 7~ 1 ai2- The derivative operators then as- 
sume the isotropic form, d^ x > — ► c?i — i02j f^ 2 -* ^82 — id%. 
(It must be kept in mind, however, that in the two sectors 
different scaling operations were performed. This type of 
scaling is only meaningful for those elements of the theory 
which do not involve a coupling between neighbouring 
nodes. We will therefore undo the scaling immediately 
after the individual sectors have been bosonized.) 

The bosonization of the free Lagrangian £0 then leads 
to the sum W[M X ] + W[M 2 ] of two WZW actions, where 
W[M n ] represents the sector (n,n) and M n (n = 1,2) 



are two independent fields taking values in the Rieman- 
nian symmetric superspace of type C\D associated with 
the supergroup OSp(4|4). From the symmetry of the 
orthosymplectic currents, J„ = ipnipn — —TiJn T i an d 
Jn = 4>n^n = — T2^T2, the field M n inherits the prop- 
erty 

A!' 1 = T 2 M S n T2 . 

As before, r k (k — 1, 2, 3) denotes the Pauli matrices act- 
ing in node space (presently, the two-component space 
underlying the definition (|2j])). To prepare the bosoniza- 
tion of the remaining operators, including those that cou- 
ple the nodes, we now undo the above scaling operation. 
This is achieved by making the replacements 

W[Mi] -» W[Mx,i] , W[M 2 ] -> W[M 2 , 7 _1 ], 

where 

W[M, 7] = ^|-r[M] - y^- J dx ld x 2 x 

x STr (^diM^diM + 1 d 2 M- 1 d 2 M) 

is the anisotropic analog of the WZW action above. No- 
tice that, owing to its topological character, the WZW 
term T[M] is not affected by the scaling operation. 

We now bosonize the remaining content of the theory, 
by making use of the dictionary. We obtain 

Coo = Y (^n T 3^n)(V'n' r 30n) «-» 
n=l,2 

<-> Y STr (t 3 J„) STr (r 3 J„) , 

71=1,2 

O'on = (^i^iH^^ + W^iXV^a^)^ 

~ STr (r 3 Ji) STr (r 3 J 2 ) + STr (r 3 Ji) STr (r 3 J 2 ) , 

°™ = 2 Y Y ('0n T fcV'™)('0^fcV'r l ) *-> 
n=l,2 fc=l,2 

~ \ J2 J2 STr ( Tfc J ") STr ( Tfc J «) - 

71=1,2 fc=l,2 

iCrO = ^2^1 + Y ■01 T fc0 ; 20'2 T fc'01 
fe=l,2,3 

<-> r 2 STr AhM 2 + r 2 Y STrM 1 T k M 2 T k , 

fe=l,2,3 

£ B = -2iE ]T " ~2ij STr M n . 

71=1,2 71=1,2 

The sum 



5 = W[M ul ] + W[M 2l 7 _1 ] + / d 2 r £ port , (26) 
£ P crt = .gCoo + 5'Coo + 9™®™ + g^oO-iio + (27) 

represents our final result for the low-energy effective ac- 
tion of the disordered d-wave superconductor: two WZW 
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actions, which are coupled by a number of marginally 
relevant operators due to inter-node and intra-node scat- 
tering of the quasi-particles. Note that the rotational 
symmetry of Euclidean space (or Lorentz invariance of 
Minkowski space) is broken by the term O' Q0 , which is a 
relict of the order parameter symmetry d x 2_ y 2 of the su- 
perconductor. In the following subsections we are going 
to discuss some physical consequences of the field theory 



A. Hard scattering, T-invariance (class CI) 

We first discuss the case where inter-node scattering is 
present, and is strong. By construction of the bosonized 
theory, the RG flow equations for its couplings coincide 
with those for the perturbed Dirac theory. From jl9| ) 
we then know that all couplings, including g^o, increase 
under renormalization. While Coo, Cooi an d Otttt are 
current-current perturbations, the operator O^o is seen 
to act as a "potential" . We expect such a perturbation 
to make some of the fields massive and remove them 
from the low-energy theory. To elucidate this effect, we 
parametrize the fields M n as follows: 



M n = exp 



A n + A r , 



B n 



-At 



A,- 



(n = l,2) 



The requirement M~ a = t 2 M^t 2 € OSp(4|4) is satisfied 
by imposing the conditions (A^) s = ±A^, B^ = B n , 
and = C n for n = 1,2. In this parametrization, the 
expansion of O,ro around unity reads 



? 2 O 7r0 = STi\(A+) 2 + (A+) 2 

+ (A- + A-) 2 + Bid + B 2 a 



D, 



It is now extremely important that, by construction [ 
of the Riemannian symmetric superspace of type C 
this potential energy is bounded from below by zero. The 
low-energy configurations of the fields Mi,M 2 are those 
that minimize O^o, which implies 
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By inserting the solution of this equation into the 
parametrization for M n , we see that the low-energy fields 
are scalar in node space: 



Mi = 



M 
M 



= Mo 



where M, a supermatrix of size (2 + 2) 2 , is subject to 
the condition M~ x = M s . The last equation defines the 
orthosymplectic supergroup OSp(2|2) with "orthogonal" 
bosons and "symplectic" fermions. (Technically speak- 
ing, M takes values in a Riemannian symmetric super- 
space of type D\C inside OSp(2|2).) Inserting the con- 
strained form of Mi and M 2 into the operators Oqq, O' 00 , 



and Otttt, we see that all of these vanish. The effective 
action thus reduces to 

S cS = 2W[M; 1 ] + 2W[M~ 1 ; 1 ~ 1 } 



7 + 7" 
8tt 



d r STr d^M~ d^M . (28) 



This is isotropic and independent of the disorder 
strength, and its ff part coincides with the replica field 
theory written down by Senthil et al. |§. The WZW 
terms have canceled because r[Mi] + T[M 2 ] = for 
Mi = M 2 -1 . 

In view of recent statements to the contrary flfi)] , we 
emphasize that the cancellation of WZW terms is a ro- 
bust feature that does not depend on any specific model 
assumptions made (as long as the disorder is generic, 
placing the model in class CI). The basic mechanism 
behind the cancellation is the dependence of the WZW 
term on parity or, equivalently, the choice of orientation 
of two-dimensional space: its sign gets reversed by the 
transformation x\ <-> x 2 . Although this dependence by 
itself does not forbid the presence of a WZW term for 
parity-invariant systems (as the sign change can be ab- 
sorbed by a target space isometry M i— > A/ -1 ), it does so 
for a d-wave superconductor in zero magnetic field. The 
crucial fact here is that the pure Dirac Hamiltonians for 
the two pairs of nodes (1, 1) correspond to opposite ori- 
entations (they map onto each other by the transforma- 
tion x\ <-> x%). Therefore, in the bosonized theory they 
are represented by WZW actions with topological cou- 
pling constants that still carry opposite signs, provided 
that uniform conventions for assigning Dirac bilinears to 
WZW fields are in force. (In the above treatment, we ar- 
ranged for the signs to be identical by choosing our con- 
ventions to be different for the two nodal sectors. This 
was done for notational convenience.) When any kind 
of scattering (consistent with the generic symmetries of 
the system) between neighbouring nodes is turned on, 
the two WZW fields Mi and M 2 become locked in the 
low-energy theory, now by Mi = M 2 , and the topological 
couplings with differing signs inevitably add up to zero. 

The model (|2g|), with the WZW term being absent, 
is called the principal chiral non-linear sigma model on 
OSp(2|2) (more precisely, on a Riemannian symmet- 
ric superspace of type D\C, which is a supermanifold 
based on the direct product of the non-compact space 
R+ = SO(2,C)/SO(2) with the compact group Sp(2)). 
Referring for a more detailed discussion of its replica ana- 
log to Ref . H , we here review only one salient feature of 
this theory. According to Friedan jnj , the one- loop beta 
function of any 2d non-linear model is determined by 
the Ricci curvature of the target manifold. For symmet- 
ric superspaces, just like for ordinary symmetric spaces, 
the Ricci curvature is proportional to the metric tensor, 
and in the present case is easily shown to be positive 
| |72[ . This means that the beta function is positive, and 
the coupling t — 871-/(7 + 7"*) therefore increases under 
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renormalization. By plausible extrapolation to strong 
coupling (t — > oo), one then expects the theory to be 
in an insulating phase with vanishing (spin) conductance 
u s = 2/(wt) — > 0. This phase was named the "spin insu- 



lator" 



Ref. 



The local density of states v of the spin insulator is pre- 
dicted to vanish linearly, v(E) ~ \E\, at ultra-low ener- 
gies, by a scaling argument due to Senthil and Fisher jl3) . 
On the basis of the supersymmetric field theory ( p8| ) , we 
can phrase the argument as follows. The growth of the 
coupling t under renormalization means that, while there 
is asymptotic freedom (accompanied only by small field 
fluctuations) at short scales, the fluctuations of M grow 
strong for large wave lengths. The scale for crossover 
from weak to strong coupling is set by the localization 
length, £ - e^" 1 §. Let us therefore partition the sys- 
tem into blocks of linear size £. To mimic the crossover 
between short-range order and long-range disorder, we 
take the fields on different blocks to be independent, and 
on each individual block to be spatially constant. By 
doing the field integral in this simple approximation we 
find 

v(E) = v f{Ev Q e) , 

where vq is the SCBA density of states, which depends 
only weakly on energy, and the scaling function f(x) 
(with asymptotic limit /(oo) = 1) has the small- a; ex- 
pansion f(x) — irx/4 + 0(x 2 ). Thus the local density 
of states goes to zero linearly at E = 0, with the char- 
acteristic energy scale being given by (i^o£ 2 ) -1 ) the level 
spacing for one localization volume. 



B. Soft scattering, T-invariance (class AIII) 

We now set the inter-node couplings 3 W7r ,3u-o to zero. 
What then remains is the free-fcrmion theory, as rep- 
resented by the sum of WZW actions W[Mi\ + W[M 2 ], 
marginally perturbed by the operator gOoo + g'0' w . Two 
observations simplify the analysis of this theory. Firstly, 
we may omit the term O' 00 when calculating the density 
of states. The physical reason is that in the absence of 
inter-node coupling we can project on one pair of nodes, 
say (1, 1), and take the other one into account by mul- 
tiplying the final answer with a factor of 2. Thus, we 
may set ip2 = "02 = and drop the field M%. Alterna- 
tively, the reduction can be seen by reasoning within the 
field-theoretical formulation, where it comes about be- 
cause the functional integral giving the contribution to 
the density of states from the pair (1,1) has an invari- 
ance under BRST transformations acting in the sector 
(2, 2). Secondly, we may rescale the coordinates to make 
the field theory isotropic. We are then facing an isotropic 
WZW model with a current-current perturbation: 



From our earlier computation of the RG beta functions 
we know that the coupling constant g does not flow. 
This can be verified directly from the WZW model by 
observing that the OPE of the Abelian current STrT3J 
with itself does not generate any UV singularities. Thus 
the above action is a fixed point for the renormalization 
group, and the theory is conformally invariant. 

What are its properties? A quick approach to the an- 
swer is to recall that the physical problem we are deal- 
ing with here is the much studied |S 73 3{| problem of 
Dirac fermions in a random Abelian vector potential. 
This problem can doubtless be analysed via the above 
action based on OSp(4|4). Unfortunately, that formu- 
lation does not seem to be the most efficient one to use 
and, in any case, it is not the one used in related previous 
work. What forced us to introduce the OSp(4|4) target 
space was the fact that we do not know of any other way 
of bosonizing the perturbations that couple the neigh- 
bouring nodes. In the present problem, however, where 
the coupling between the nodes is absent, the symmetry 
group consists of two copies of GL(2|2), the subgroup of 
elements h of OSp(4|4) with the property hr^h^ 1 = T3. 
The excess of field degrees of freedom in OSp(4|4), as 
compared to GL(2|2), means that there is some redun- 
dancy. An efficient method of solution will try to avoid 
this redundancy, and can avoid it as follows. 

In standard non-Abelian bosonization (without boso- 
nic ghosts), two schemes are distinguished: bosoniza- 
tion takes the free-fermion theory into a level-one WZW 
model either over the orthogonal group 0(2n), or over 
the unitary group U(n). While the two bosonization 
schemes are equivalent in the free limit, the 0(2n) version 
has the advantage of exposing the full set of symmetries 
of the free theory of n Dirac fermions (or 2n Majorana 
fermions), whereas the U(n) version is better adapted 
to the treatment of certain perturbations. The latter 
scheme was developed particularly by Affleck [(74) in his 
treatment of the Id Hubbard and Heisenberg models. 

By straightforward transcription of Ref. [Q, we can 
generalize the U(n) bosonization scheme to include a 
bosonic ghost system. The resulting theory is again a 
supersymmetric WZW model, with an action functional 
of the same form as before, but the field now takes val- 
ues in the supergroup GL(2|2) (or, rather, in a Rieman- 
nian symmetric superspace of type A\A inside GL(2|2)). 
Deferring the details of this generalization to Appendix 
[B|, we here take a short cut. To arrive at the bosonized 
theory in the desired form, we start from the action func- 
tional ( |29| ) and simply insert for M\ the reduced expres- 
sion 



I M 



where M £ GL(2|2) is parametrized, for example, by 



W[M]+g / d 2 rSTr(r 3 J)STr(T 3 J) 



(29) 
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The action functional for M then becomes 

iT[M] 
12tt 



S[M] =--!- / cPrSTtduM- l daM- 
8n J 



4- / d 2 rSTr(M~ 1 dM) STr(M<9M _1 ) . (30) 



The same result is obtained by direct application of the 
GL(2|2) bosonization rules to the supersymmetric Dirac 
theory with the marginal perturbation g J OaocPr. The 
numerical (or bosonic) part of this action has the impor- 
tant property of being bounded from below (for g not in 
excess of a certain critical value), and hence leads to a 
well-defined functional integral. (Positivity of the action 
is ensured by the construction of the field manifold as a 
Riemannian symmetric superspace.) 

Past solutions of the physical problem of Dirac 
fcrmions in a random Abelian vector potential have used 
the equivalent representation by a replica sine-Gordon 
model (H), or diagonalization of the quantum Hamilto- 
nian by a Bogoliubov transformation J73[ , or the study 
of Kac-Moody current algebras with U(ljl) x U(l|l) sym- 
metry . A closely related model has been investigated 
via sigma model techniques in j75|. Here we are going 
to take a different approach: we will solve the problem 
by direct manipulation of the functional integral. Fol- 
lowing Knizhnik and Zamolodchikov |76| , this is possible 
by exploiting the local (infinite-dimensional) symmetry 
G L x G R (with G = GL(2|2)). For better readability the 
technical details of this computation have been relegated 
to Appendix [b|, and we now just give a summary of the 
most important points. 

The basic idea is to study the response of the functional 
integral to field variations of the form 

S x M(r) = -X(r)M(r) , 

where X(r) takes values in the Lie algebra gl(2|2). As 
usual, variation of the action yields the Noether current: 



6 X S = - I d 2 rJ gx , 



Jx = STriXMdAr 1 ) + — STr(AT)STr(M9M" 1 ) . 

7T 

From the stationarity condition 6x S = 0, one infers that 
J is a holomorphic conserved current, i.e. BJa — for 
spatially constant A and on solutions of the equations of 
motion. In addition, the theory has a conserved current J 
which is antiholomorphic: d J a — 0. The expression for it 
is determined by the change of S in response to variations 
5yM(r) = M(r)Y (r). In the following we concentrate on 
the holomorphic sector. 

By exploiting the invariance of the functional expec- 
tation value (...Jb(0)) (where the ellipses indicate ad- 
ditional operator insertions) w.r.t. the variable change 
M -> c~ x M = M + S X M + ... with S X M = —AM, one 
deduces the operator product expansion (OPE) for the 
currents: 



M , um .f!*m + $*-M + m , (31 , 

z z z 
f(A, B) = -STrfAB) + — STr(A)STr(B) . 

Like any operator product expansion, it is to be under- 
stood as an identity that holds under the functional in- 
tegral sign. The dots indicate terms which remain finite 
in the limit z — > 0. 

The behaviour of correlation functions under scale 
changes or, more generally, under conformal transforma- 
tions z 1— > s(z), is determined by the OPE of the fields 
with the holomorphic component of the stress-energy- 
momentum tensor, T{z). To obtain the latter, one de- 
mands that the OPE between T(z) and the currents 
starts as T(z)Ja(0) — Ja(0)/z 2 + which expresses the 
fact that Ja{z) has conformal dimension (1,0). Assum- 
ing T(z) to be of the Sugawara form und using Eq. (|3l]), 
one then finds 

T(z) = —J^J^z) + JUL J e (z)J e {z) . (32) 

The notation here means the following: choosing some 
basis {e^} of the Lie superalgebra gl(2, 2), we set Jj = J ei , 
and Kij = — STr etej, and we raise the indices of the met- 
ric tensor by K^Kjk = 5 l k . The current J e , whose square 
gives the second term in the formula for T (z), represents 
the generator e = id (unit matrix). 

The observable we are interested in is the local den- 
sity of quasi-particle states at low energy. Recall from 
(|15| ) that a finite energy E is accounted for in the field 
theory by a perturbation -HE /d 2 r(*f* 1 + ^f^j). 
This bosonizes to -2i(E/£) J d 2 r STr(M + M^ 1 ). More- 
over, the local density of states at r = is given 
by the expectation value of the bosonic (or fermionic) 
part of ^(O)*!^) + *f(0)*j(0), which bosonizes to 
(2^)- 1 Tr (M(0) + Af(0) _1 ). 

Thus, the density of states is determined by the ex- 
pectation value of the fundamental field M and its in- 
verse M _1 . To work out its scaling behaviour, we need 
the conformal dimension of M , which is the number Am 
multiplying the leading singularity in the OPE of T{z) 
with M(0): 

T(z)M{0) = A M M(0)/z 2 + ... . 

Given the Sugawara form ( |32| ) of the stress-energy- 
momentum tensor, the value of Am follows from the OPE 

J A (z)M(Q) = -AM(0)/z + ... , 

which reflects the fact that, by construction, J a repre- 
sents the generator of left translations M 1— > e~ A M. Us- 
ing that the quadratic Casimir in the fundamental rep- 
resentation of GL(2|2) vanishes (tfieiej =0), one finds 
Am = 5 — g/it- The total scaling dimension of M is the 
sum (from holomorphic and antiholomorphic sectors) 
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A = A M + A M = 1 - 2g/ir . 
To calculate the local density of states, v{E), we add 

-HE I d 2 rSTr(M + M" 1 ) + ATr (M(0) + M(O)" 1 ) 



to the action functional and differentiate with respect 
to the source A at A = 0. If the short-distance cutoff is 
raised by a renormalization group transformation r t— > 6r, 
the energy scales as E i— > b 2 ~ A E and the local source as 
A h- ► fe~ A A. Consequently, v{E) obeys the homogeneity 
relation 

v{E) = b- A is{b 2 - A E) . 
The solution is a power law 

Thus the density of states is expected to vary alge- 
braically with a non-universal exponent that depends on 
the disorder strength g. This is in agreement with a re- 
sult first found by Nersesyan, Tsvelik, and Wenger ||. 

Let us mention in passing that the critical line for Dirac 
fermions in a random Abelian vector potential has been 
argued to be unstable with respect to an infinite 
number of relevant perturbations. In our opinion, this 
instability is not a true effect but disappears when the 
theory is formulated on a tar get manifold with Rieman- 
nian structure (see Appendix B 1). 



C. Realistic scenario? 

Having reviewed the two extreme cases of strongly cou- 
pled and completely decoupled nodes, it is pertinent to 
address what kind of scenario will prevail under realistic 
conditions. By the very nature of the question, no clear- 
cut answer exists. Yet assuming that the basic modeling 
by the microscopic Hamiltonian (||) has some finite over- 
lap with experimental reality, we believe that two basic 
scenarios are conceivable: 

Systems with a significant amount of large-momentum 
transfer scattering on the microscopic level. Here the 
adjective "significant" refers to a regime in which the 
strength of the inter-node coupling is in excess of the 
other characteristic energy scales of the problem (such as 
temperature, experimental resolution, etc., which couple 
to the action via the energy parameter E) . In that case, 
relative fluctuations of the two fields Mi and Mi are neg- 
ligible, and it is only the isotropic content of the theory, 
described by the action (p8|), that matters. 

A more complicated situation arises in systems whose 
microscopic Hamiltonian contains only soft scattering: in 
systems where the inter-node scattering is a small pertur- 
bation, we run into a crossover scenario. The bare theory 



is given by two WZW actions, weakly perturbed by inter- 
node scattering. As we saw, renormalization makes the 
inter-node couplings grow. Therefore, the theory flows 
to strong coupling on asymptotically large scales and the 
fields for the different nodes will again be locked. How- 
ever, that fixed-point regime may be preceded by a large 
intermediate region in which the critical properties of the 
WZW theory prevail. In that case we expect for the DoS 
a power law which is both non-trivial and non-universal, 
as predicted by (|33|). 

Summarizing, the d-wave system renormalizes to a uni- 
versal limit described by the rotationally and parity in- 
variant action (p8[). There may, however, be an extended 
crossover region in which the system of decoupled single- 
node actions dominates. The question of which type of 
scenario prevails in real and numerical experiments is sys- 
tem specific and not for the present approach to decide. 



D. Soft scattering, broken T-invariance (class A) 



We have argued in Section II B that perturbations 
breaking time-reversal invariance reduce the symmetry 
of the (i-wave superconductor without inter-node scat- 
tering from class A lll t o A. In the mixed state, as we 
recall from Section II A , T-breaking is due primarily to 
the supercurrent flow, which gives rise to a scalar po- 
tential in the quasi-particle Hamiltonian projected on a 
single node. This issue is resumed in the present subsec- 
tion, where we address the perturbation of the class AUI 
Dirac Lagrangian (for the pair of nodes 1,1) by 

3 is the second component of the su- 
percurrent velocity. According to the bosonization rules 
derived in Appendix 0, this perturbation transforms into 



where v s = 



l STrof ■ M 



As before, M takes values in (the Riemannian symmet- 
ric superspace of type A\A inside) GL(2|2). Thus the 
complete action functional is 



M~ 



where S[M] was defined in Eq. (|30|). The perturbation 
has positive scaling dimension 2 — A = f + 2g/7r and is 
therefore relevant. 

We expect such a perturbation to cause a reduction 
of the field manifold. Unlike the symmetry-breaking 
terms we considered earlier, the present one is imaginary- 
valued, giving rise to oscillations in the functional inte- 
gral. The argument for field reduction therefore can- 
not be the rigorous kind of potential-energy argument 
that applied to the crossover from class Alll to CI. 
What we must rely on now is stationary phase, i.e. the 
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fact that contributions from field configurations with in- 
dependent phases tend to cancel out. Starting from the 
identity element M = 1 and conjugating by a symme- 
try transformation for class A, M i— > TM<j^T~ 1 g^° 
with T 6 GL(2|2), we get a stationary-phase manifold 
parametrized by 

Mag 3 = V.; M 1 = TafT 1 = Q . 

By construction, the above T-breaking perturbation van- 
ishes on this manifold, 

STrof^M + M- 1 ) -> . 

As for the other terms in the Lagrangian, it is straight- 
forward to show that they reduce to 

STrdpM-%M STr d^Qd^Q , 
STr(M<9M _1 ) STr(M _1 3M) -> . 

The fate of the WZW term T[M] under the field reduc- 
tion is less trivial. Doing the same steps as in Section 7 
of Ref. Q], we find that T[M] reduces to a topological 
theta term on the manifold of Q-matrices: 

LZn J 

Cto P = -^STrQd tt Qd u Q, 
lorn 

with 8 = ±7r. The ambiguity in the sign of 8 results 
from the multivalucdness of the WZW term: to calculate 
T[M] one must extend the field M to a 3-space bounded 
by the two-dimensional position space, and there is no 
unique way of doing that. However, if the position space 
is closed (as is implied by its being a boundary) the am- 
biguity in the sign of 8 has no consequence. In that case 
it can be shown that the topological term is a winding 
number or, in other words, the integral of the topologi- 
cal density £ t op over any closed surface takes quantized 
values in 2niZ, so that a shift of 8/2ir by any integer n 
is unobservable in the functional integral. For the pair of 
nodes (2, 2) the situation is similar. 

To summarize, we find that the low-energy effective 
action for the pair of nodes (n, n) of a c£-wave supercon- 
ductor in symmetry class A is 

S^[Q]=±S top [Q] 

--^ J d 2 r STr (j^d^diQ + j s -d 2 Qd 2 Q) , 

where s\ = 1 and S2 = — 1. In the isotropic limit 7 = 1, 
this action functional is familiar from the integer quan- 
tum Hall effect, where it is known as Pruisken's action 
Ifrst with longitudinal electrical conductivity o~ xx = 1/tt 
and Hall electrical conductivity a xy — 8/2-k = ±1/2. 
The physics associated with this action is (i) a finite and 
smooth density of states at E = 0, and (ii) critical be- 
haviour in the transport coefficients. 



At first sight, one might be suspicious of our nonzero 
result for the coupling a xy = 8/2ir. What is puzzling 
about it is this. Since the Hall conductivity transforms as 
a pseudoscalar, i.e. has its sign reversed by the exchange 
of Cartesian coordinates x y, a nonvanishing Hall re- 
sponse breaks parity. On the other hand, the massless 
Dirac theory is invariant under parity (if the transfor- 
mation x <-> y is accompanied by the exchange of the 
left-moving and right-moving fields) and so is the WZW 
model (if x <-> y is combined with M <-» M -1 ). It goes 
without saying that a parity-invariant theory cannot gen- 
erate pseudoscalar observables. On these grounds, recent 
papers have stated J79| , ^6|j80| that the thermal Hall con- 
ductivity of a d-wave superconductor in the mixed phase 
vanishes in the linear (or Dirac) approximation. 

Does it follow from this symmetry argument that 
the reduction to the nonlinear sigma model and/or our 
nonzero result for the coupling a xy must be false? It does 
not. Indeed, recall the invariance of expid J d 2 r£ top un- 
der shifts 8 — ► 8 ± 2ir, when the position space has no 
boundary. This invariance means that a xy = 1/2 and 
Cxy — — 1/2 give equivalent field theories for the case of a 
closed system. Thus, although a xy is a pseudoscalar, the 
nonlinear sigma model at a xy — ±1/2 enjoys the prop- 
erty of being parity- invariant, and there does not exist 
any conflict with the parity invariance of the massless 
Dirac theory (or the WZW model). 

At the same time, it would be a preposterous proposi- 
tion to claim that Pruisken's nonlinear sigma model for 
the integer quantum Hall effect (IQHE) at a xy = ±1/2 
describes a parity-invariant physical system with vanish- 
ing Hall response! How do we then reconcile parity in- 
variance at a xy — ±1/2 with IQHE phenomenology? The 
discrepancy is resolved by observing that the Hall con- 
ductivity is not a Fermi-edge quantity, which is to say 
that it is not determined by the Green's functions at a 
single energy E = Ep only, but is given by a sum over 
energies. This feature of the Hall conductivity, which 
distinguishes it from the longitudinal conductivity a xx , 
poses a difficulty to the nonlinear sigma model formula- 
tion: since the field theory is derived for a fixed energy, all 
the information needed to reconstruct the Hall response 
in the bulk cannot be contained in it. In other words, the 
topological coupling a xy of the nonlinear sigma model 
does not coincide with the total Hall conductivity of the 
bulk, and hence the parity invariance of the field theory 
at a xy — ±1/2 does not imply that the Hall response of 
the quantum Hall system vanishes. 

Nonetheless, as was pointed out by Pruisken Jr8[ , it is 
possible to monitor the Hall response field-theoretically 
by studying systems that have a boundary. There ex- 
ist chiral boundary currents at the edges of quantum 
Hall systems, and the topological term of the nonlinear 
sigma model (for a fixed energy E = Ep) does cap- 
ture the physics of these boundary currents. In fact, 
when a boundary is present, the invariance under shifts 
o~xy o~xy ± 1 no longer holds, parity symmetry is bro- 
ken for all nonzero values of o~ xy including a xy = ±1/2, 



2G 



and the topological term gives rise to long-range current- 
current correlations at the edges of the system. The mes- 
sage from this is that, in order to get an unambiguous 
field-theoretic view of thermal and spin Hall transport in 
a d-wave superconductor, we are well advised to consider 
geometries with a boundary. 

Having understood that, we run into another difficulty: 
A consistent extension of a single WZW model to a sys- 
tem with boundaries is not possible! Instead, a geometric 
construction detailed in Appendix ^] will show the follow- 
ing: 

• The modeling of a system S with boundary <9£ re- 
quires, at least, two WZW fields Mi and M 2 . 

• These fields are coupled at the boundary through 
a condition (Afi) _1 | as = AfaLj, (plus additional 
conditions on the derivatives of the fields to ensure 
current conservation.) This mechanism reflects the 
fact that specular reflection at the boundaries of a 
finite systems couples quasi-particle species of dif- 
ferent parity. (Specular reflection inverts one of the 
momentum components while leaving the other in- 
variant: a parity operation.) 

• As far as the Hall conductivity of the coupled 
system is concerned, everything depends on the 
masses mi,2 of the two field species: (i) if both 
fields are massless, they enter the theory with 
equal weight and the Hall conductivity vanishes 
by symmetry, (ii) In highly anisotropic cases, e.g. 
mi S> m 2 , the field M\ decouples from the low en- 
ergy sector of the theory. Under these conditions, 
a xy = ±1/2 depending on which field has become 
heavy, and the sign of rni.2- Finally, (iii) if the 
masses \m\\ ~ |m2| are approximately equal, it is 
the signs of mi <2 that determine the Hall conduc- 
tivity. If the signs are the same, we expect to get 
<T xy = ±1; if they differ, we expect a xy = 0. 

What can we infer from these observations about the 
problem at hand? Recall that the two pairs of conju- 
gate nodes (1, 1) and (2, 2) are represented by the WZW 
fields Mi and M2 respectively, with our conventions be- 
ing such that the total WZW term is T[Mi\ + T[M 2 ]. 
Given these low-energy effective fields, we impose the 
boundary condition (Mi ) ~ 1 1 as = M 2 1 as . (Alternatively, 
the boundary conditions could mix in high-energy fields 
corresponding to excited quasi-particle bands of the su- 
perconductor.) From what we said before, this boundary 
condition is forced by the requirement that the WZW 
functional be mathematically well-defined. If the parity 
invariance of the microscopic Hamiltonian is broken by 
the presence of a magnetic field in the form of vortices 
and by the supercurrent flow, there exists no fundamen- 
tal symmetry principle that would protect the fields from 
acquiring (small) masses mi_ 2 in the low-energy effective 
theory. 

More detailed analysis would then be required to de- 
cide which scenario (i-iii) applies in the case of a realistic 



d-wave superconductor (for the most complete investi- 
gation to date, see @). However, as far as theory is 
concerned, let us re-iterate that the WZW model and 
the original Dirac fermion model represent equivalent 
descriptions of the system. In view of the consistency 
problems outlined above, and their resolution by an ex- 
tended theory, one expects that a faithful description of 
Hall transport in the Dirac fermion language, too, has 
to involve all four nodal fermion species simultaneously. 
Therefore, symmetry-based statements that the Hall con- 
ductivity of a single Dirac fermion system vanishes p(| 
should be taken with a grain of salt. 

To avoid confusion, we elaborate on one other subtle 
point that may seem paradoxical. Recall that what we 
set out to calculate was the single-particle Green func- 
tion, giving the density of states. Since the field theory 
we have derived for it contains massless (Goldstone) fields 
Q, one might be misled into thinking that our theory pre- 
dicts singular corrections to the density of states. In the 
present instance, this is not so. The reason why we ob- 
tained massless fields is that, in the course of the deriva- 
tion, we chose to carry out a doubling of field variables 
(introducing charge conjugation space, so as to accom- 
modate the symplectic symmetry of the Gorkov Hamil- 
tonian, which becomes important as soon as the coupling 
between nodes is turned on). This means that we are 
effectively working in the two-particle sector. Put dif- 
ferently, the present theory involves all of the fields for 
the nodes (1,1), although the basic quasi-particle Hamil- 
tonian connects only half of them. As a result, when 
we compute the density of states the source fields cou- 
ple only to a subset of the Grassmann fields. Thus there 
is a global fcrmionic symmetry, which causes the field 
integral to collapse by the usual BRST mechanism and 
yield a vanishing correction to the density of states. The 
field theory stays alive only for the two-particle Green 
function, viz. transport. 

E. Hard scattering, broken T-invariance (class C) 

The final case to consider is the simultaneous presence 
of inter-node scattering and breaking of time-reversal in- 
variance. Having switched to the GL(2|2) bosonization 
scheme for convenience of presentation and solution in 
the previous two subsections, we must now return to the 
more general OSp(4|4) scheme in order to accommodate 
the inter-node couplings. The T-breaking perturbation 
by a supercurrent flow is rewritten in terms of ip n , tp n as 

i^Vf ® r 3 )vK = iSTr (a 3 cc ® r 3 ) + . 

By the OSp(4|4) bosonization rules, this bosonizes to 

U~ 1 STr (af' <gi r 3 ) (M 1 + M 2 1 ) . (34) 

Now recall that the effect of inter-node scattering is to 
lock the fields (Mi — M^ 1 ) and to reduce Mi to a scalar 
in node space, 
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with M s = M . Inserting this expression for Mi and 
using Tr r 3 = 0, we find that the perturbation disappears. 
Thus, in the linear approximation used here, the break- 
ing of time-reversal symmetry by the supercurrent flow 
ceases to be effective in the node-locked limit. 

While this may seem surprising at first, it is not hard 
to comprehend in physical terms. The perturbation 
we are considering represents the Doppler shift due to 
the motion of quasi-particles relative to the supercur- 
rent flow v s . Because the node 1(2) and its conjugate 
1(2) are located at opposite momenta, the corresponding 
quasi-particle states experience opposite Doppler shifts. 
Therefore, scattering between opposite notes causes the 
Doppler effect to alternate randomly in time. It is clear 
that the node-locked limit will be attained over scales 
large enough in order for a quasi-particle to undergo 
many scattering events between the nodes. On such 
scales, and in modes of long-range quantum interference 
that are stable with respect to disorder averaging, the 
Doppler shift averages to zero. 

However, this does not mean that the presence of the 
supercurrent flow is of no consequence. At the free- 
fermion point, the field M has scaling dimension one, 
which makes the perturbation (Q) more relevant in the 
RG sense than the marginal perturbations that represent 
inter-node scattering. (Precisely speaking, this is true 
only for wave lengths smaller than the spacing between 
vortices. On large scales, the spatial oscillations of the 
supercurrent velocity cause the perturbation to average 
to zero.) We should therefore evaluate the consequences 
of the perturbation ( |34j) before locking the nodes. Doing 
so, and appealing to the same argument as in the previous 
subsection, we expect the perturbation ( |j) to reduce the 
field M n (n = 1,2) to stationary-phase manifolds which 
now are given by 



Mi 
M 2 



00 



00 
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(^ c ®T3)r 2 (a3 CC ®r 3 )r 2 - 1 



with T n 6 OSp(4|4), the diagonal subgroup of the chiral 
symmetry group of the C\D WZW model. 

Next, we take the effects of inter-node scattering into 
account (with the plausible assumption that reduction of 
the field manifold does not change the relevant nature 
of the perturbations from inter-n ode s cattering.) By the 
mechanism described in Section IV A , the fields will be 
locked: 



T 2 = 



T 
T 



where T now belongs to the other variant of OSp(2|2) 
(orthogonal bosons, symplectic fermions). By inserting 
the locked fields into the WZW action, we arrive at a 
non-linear sigma model with action 



S[Q] = -\ J d 2 rSTrd^Qd^Q 
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t=— — j, (35) 
7 + 7 1 

where Q = TafT parametrizes a Riemannian sym- 
metric superspace of type D\C. Note that in the lin- 
ear order of approximation used, there is no topological 
theta term in this non-linear sigma model. Precisely the 
same field theory (without a definite value for the cou- 
pling t) was obtained in |67| from a more phenomeno- 
logical approach. The known physics of this system is a 
Zeeman-field dependent weak-localization correction to 
the thermal conductivity |67j], and a density of states 
that vanishes as E 2 at ultra-low energies |0 . 

Rather than elaborating on these properties, we turn to 
another way of breaking time-reversal symmetry, which 
is to make the order parameter of the superconductor 
complex by including an imaginary component. Moti- 
vated by recent developments that were reviewed in the 
introduction, we consider an id xy component with mo- 
mentum dependence iXsin(k x a) sm(k y a). In the Gorkov 
Hamiltonian, such a term enters as the Pauli matrix oi 
acting in particle-hole space. If we again rotate by e l7Tcri ^ 4 
and linearize around the nodes, the term becomes pro- 
portional to (73 and acts as a Dirac mass term, with the 
mass parameter being positive for one pair of nodes, say 
(1,1), and negative for the other pair. In the field the- 
ory, upon bosonization, the perturbation is represented 
by the operator 

' A ' (36) 



■STr of 3 (M x 



M2 1 ) 



which is relevant at the free-fermion point. The observa- 
tion we made earlier for the case of a supercurrent flow 
applies here as well: if we assume the limit of locked 
fields (Mi = M^ 1 ), the perturbation vanishes. This is 
again plausible, since the Dirac masses resulting from 
sin(fcj;a) sin(fc y a) carry opposite signs for neighbouring 
nodes, which, as we recall, transform into each other by 



k~r or kn 



-k y . Hence the Dirac mass of a 



quasi-particle subject to inter-node scattering is zero on 
average. As before, however, the correct procedure is first 
to deal with the relevant mass perturbation, and only 
then to impose field locking. Proceeding in this way, we 
would again be led to the action ([35]) for the low-energy 
effective theory. 

We should, however, be suspicious of the last result, as 
it does not conform with what is expected on symmetry 
grounds. The point is that the action of the non-linear 
sigma model (35) is invariant under a parity transfor- 
mation x\ 1 ► x\ and xi > —X2, whereas the perturbed 
WZW theory we are starting from is not. Indeed, if we 
combine the parity transformation with field inversion 
Mj 1— > M" 1 , so as to preserve the WZW action, the sign 
of the perturbation gets reversed: 



STr £7<f°Mf 1 = STra^Ml = -STr Mi a 



CO 
3 



where the last equality follows from STr(AB) = 
STr(AB)' = STr(B* J 4*) and the fact that of is odd un- 
der the orthosymplectic transpose defined by (aftpy = 
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-0*(<73 C )*- Thus, parity is broken in the initial theory, and 
it is reasonable to expect that it will still be broken in 
the final theory. 

As was first pointed out by Pruisken et al. [ [77|j78| 
in the related context of the integer quantum Hall ef- 
fect, the way to break parity invariance in the non- linear 
sigma model is to add a theta term with topological an- 
gle 9 ^ itZ to the action. From what has been said, 
we expect such a term in the low-energy effective ac- 
tion when a secondary id xy component of the order pa- 
rameter is present. To compute the angle we aban- 
don non-Abelian bosonization and revert to the stan- 
dard method: Hubbard-Stratonovich transformation and 
saddle-point approximation, followed by a gradient ex- 
pansion. (Note that, although non-Abelian bosonization 
is perfectly suited to describe the flow of the field theory 
away from the free-fermion point, it is a less reliable tool 
for understanding how the theory arrives at and flows 
into the non-linear sigma model.) The last step is read- 
ily performed using the meth od of heat kernel regular- 
ization reviewed in Appendi x [A 4 As a result we obtain 
the action presented in Eq. (Ill), with couplings 
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A summary of the physics of this model, and references 
to the relevant literature, were given in the introduction. 



V. COMPARISON WITH OTHER APPROACHES 

In this section, we comment on the second contro- 
versy mentioned in the introduction: the incompatibil- 
ity of field theory and self-consistent approaches. Sev- 
eral arguments have been put forward to explain the ex- 
istence of conflicting results for observables as basic as 
the mean DoS (finite at zero energy in the self-consistent 
approaches, vanishing in field theory). It has, for ex- 
ample, been argued that the vanishing of the DoS in 
field theory might be due to an ill-defined continuum 
limit p^ ]. On the other hand, it has been pointed out 
|S3f that the notorious appearance of logarithmic diver- 
gences in the diagrammatic results hints at an insufficient 
UV-regularization. However, in hindsight it looks like 
these, UV related aspects of the problem carry minor if 
any significance. The origin of the discrepancy between 
field theory and early self diagrammatic aproaches to the 
problem lies somewhere else: 

First of all, it is important to realize that the vanishing 
of the DoS obtained in the field-theoretical description is 
due to a mechanism which can simply not be included 
into diagrammatic approaches. Indeed, we have seen in 
Section [II that the breaking of the chiral symmetry of 
the Hamiltonian on the saddle-point level leads to the 
appearance of Goldstone modes. In early perturbative 
approaches to the problem [jl] |) , these modes had not 
been included at all. In the field theory, however, they 



are responsible for much of the physical behaviour of the 
system. In particular, it is the strong fluctuation of these 
modes in low-dimensional systems that leads to a vanish- 
ing of the DoS at E = 0. 

It is tempting to explain this vanishing by the Mermin- 
Wagner-Coleman theorem as was done by NTW. Indeed, 
a non-vanishing DoS would have meant that a contin- 
uous symmetry of the field theory was spontaneously 
broken. (As with normal disordered systems, the DoS 
plays the role of an order parameter in the field the- 
ory.) The Mermin- Wagner- Coleman theorem essentially 
states that, as a consequence of unbounded Goldstone 
mode fluctuations, spontaneous breaking of a continuous 
symmetry does not occur in dimensions d < 2. There 
is, however, a caveat with that argument, which is that 
the Mermin- Wagner-Coleman theorem explicitly refers to 
compact field manifolds. In the non-perturbative super- 
symmetry approach at hand, we are not dealing with such 
types of field manifold. This is not an academic point, as 
is shown by counterexamples: there exist cases, namely 
two-dimensional systems in class D or Dili, where the 
naive application of the Mermin- Wagner-Coleman theo- 
rem leads to the erroneous prediction of a vanishing DoS. 
(The non-linear sigma model for weakly disordered 2d 
systems in these classes predicts a zero-energy DoS which 
diverges in the thermodynamic limit J47U49].) However, 
for the systems in the classes AIII, CI, and C studied in 
the present paper, the inclusion of the Goldstone modes 
indeed leads to a vanishing DoS. 

It is very instructive to explore the phenomenon in the 
zero-dimensional case (i.e. the case of ergodic systems, 
where spatial fluctuations of the Goldstone modes are 
frozen out.) The reason why d = is a good case to study 
is that the zero-dimensional non-linear sigma model in- 
tegral can be performed explicitly (c.f., e.g., Ref. Q). 
It turns out that the results, including the vanishing of 
the DoS at zero energy, agree with the phenomenologi- 
cal random-matrix theory approach to systems of class 
C, CI, and Alll Moreover, the simplicity 

of the 0<i-case makes it particularly straightforward to 
deduce what is missing in the diagrammatic approach. 
After all, it cannot be that perturbative diagrammatic 
approaches are completely oblivious to the existence of 
the Goldstone modes mentioned above. In fact, it has 
been shown in Ref. |34j that relevant classes of diagrams 
are missed within the SCBA approach to computing the 
Green function. "Relevant" here means that the dia- 
grams in question diverge as the energy approaches zero. 
This behaviour is indicative of the fact that these dia- 
gram classes represent the perturbative implementation 
of the Goldstone modes discussed above (very much like 
the standard diffuson mode is the first-order perturba- 
tive contribution to the Goldstone modes induced by 
the spontaneous breaking of the symmetry between re- 
tarded and advanced Green functions for disordered met- 
als). More recent diagrammatic approaches [ |S6| have in- 
cluded these modes into the perturbative theory of the 
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disordered d-wave superconductor. In intermediate en- 
ergy regimes, where the IR singularity of the Goldstone 
modes has not yet become virulent, the extended dia- 
grammatic formulation represents an alternative to the 
field-theoretical approach. 

Summarizing, we find that (i) the inclusion of Gold- 
stone modes is crucial for a correct description of the 
quasi-particles of a d-wave superconductor (in particular 
its DoS), (ii) that within diagrammatics these modes are 
represented by singular diagram classes which (iii) can 
be brought under control in intermediate, but not in the 
lowest energy regimes. 

Finally we comment on work where a non-vanishing 
zero-energy DoS has been obtained in a non-perturbative 
manner. In particular, Lee Q considered the DoS on the 
background of an ensemble of impurities at the unitarity 
limit. The strength of these impurities makes a compar- 
ison with our analysis difficult. (The construction of the 
field theory crucially relies on the existence of a para- 
metric separation between the energetic extension of the 
nodal region and the much smaller width of the disorder 
distribution.) 

Another type of non-perturbative approach has been 
put forward by Ziegler, Hettler and Hirschfeld 0. In 
their work, a non-vanishing DoS was obtained on the 
basis of rigorous estimates applied to the single-particle 
Green function. Upon close inspection, however, their 
line of reasoning can be dismissed by its lack of relevance 
for superconductors. Indeed, the Hamiltonian they con- 
sider is extremely special and non-generic in class CI. It 
is formulated on a lattice and, in the notation of Ref . j| , 
reads 

where the kinetic energy —V 2 and the non-local d-wave 
order parameter A d are taken to act by 

(V 2 V)(r) = V(r + 2e x ) + ^(r - 2e x ) 
+ ^(r + 2e 2 )+#r-2e 2 ) , 

(A d V)(r) = A W>(r + ei) + ^(r - ei) 
- tp(r + e 2 ) - ip(r - e 2 )) , 

and n is a random chemical potential. The choice for V 2 
has the artificial feature that it allows hopping only be- 
tween third-nearest neighbours, which leads to a conser- 
vation law alien to superconductors: H commutes with 
the operator D0-3, where D r r ' = (— I) xi+X2 5 r>r > multi- 
plies by plus one on the sites of an A sublattice (x\ + x% 
even) and by minus one on the sites of the B sublattice 
(xi + X2 odd). Since Da 3 has two eigenvalues ±1, the 
Hilbert space decomposes into two sectors not coupled 
by H . The first one consists of all particle states on the 
A sublattice and hole states on the B sublattice, while 
the second sector contains all the complementary states. 
Without loss of generality we may restrict the Hamilto- 
nian to one sector. The particle-hole degree of freedom 
then becomes redundant — we can make a particle-hole 



transformation on the B sites, so that the first sector be- 
comes all particles and the second sector all holes — and 
the Hamiltonian reduces to 



H 



-(V 2 + m) 



A 



+ (V 2 + fi) 



A 



d 



B 



where —(V 2 + ju) L denotes the restriction of —(V 2 + fx) 
to the A sublattice. H is a generalized discrete Lapla- 
cian augmented by a random potential. It belongs to the 
symmetry class Al, the Wigner-Dyson class with (3 = 1, 
which is to say that the Hamiltonian matrix is real sym- 
metric, and the large-scale behaviour of its two-particle 
Green's function is controlled by the cooperon and diffu- 
son modes well known from the theory of disordered met- 
als. There exist no quantum interference modes affect- 
ing the single-particle Green's function. Thus the quasi- 
particle Hamiltonian of Ziegler, Hettler and Hirschfeld 
models a time-reversal invariant normal metal rather 
than a superconductor! As one would expect, it can be 
proved [|J that such a Hamiltonian has a non-zero den- 
sity of states at E = for a wide class of distributions 
of the random potential \i. However, this result tells us 
nothing about a disordered d-wave superconductor. The 
built-in Do~3 conservation law eliminates all the modes of 
quantum interference that are characteristic of the super- 
conductor and act to suppress the density of states at zero 
energy. As a corollary, we conclude that the ZHH model 
provides no test of the accuracy of the self-consistent T- 
matrix approximation widely used for superconductors. 

To preclude any misunderstanding, let us emphasize 
that the ZHH model is built on a superfluid conden- 
sate, which may of course exhibit the Meissner effect 
and other phenomena associated with superconductivity. 
However, the issue at hand is not the nature of charge 
transport by the superconducting condensate. As was 
emphasized in the introduction, we are not studying the 
condensate, but rather its quasi-particle excitations, their 
spectral statistics and their transport properties, which 
can be probed experimentally via spin or thermal trans- 
port. From this perspective the ZHH Hamiltonian, al- 
beit built on a superconducting ground state, models a 
normal metal. (More precisely, it gives the behaviour a 
thermal insulator, since quantum interference effects ul- 
timately drive the model to strong localization in two 
dimensions.) 



VI. NUMERICAL ANALYSES OF THE 
QUASI-PARTICLE SPECTRUM 

Much of the early work on quasi-particles in disor- 
dered d-wave superconductors was analytical. More re- 
cently, a number of numerical investigations exploring 
the effects of disorder scattering appeared (see, e.g., 
Refs. |^,^l|,^5|). Taking the soft and hard scattering 
regimes in turn, the present section reviews elements of 
these works, and relates them to the results disussed in 
previous sections. 
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A. Hard Scattering 

A comprehensive analysis of quasi-particle spectra 
in time-reversal invariant <i-wave superconductors with 
point-like scatterers (symmetry class CI) appeared in 
Refs. pc|] . Going beyond the mere diagonalization of 
the lattice Hamiltonian (|J) , these papers determined the 
order parameter self-consistently. Moreover, the role of a 
nesting symmetry of particle-hole type, which is present 
in the case of a half-filled band and refers to momen- 
tum transfers q = (±7r/a, ±7r/a), was explored. Without 
going into quantitative detail, the main results of these 
papers can be summarized as follows: 

• The self-consistent T-matrix approximation fails to 
correctly describe the DoS below a certain energy 
scale, the value of which increases with disorder. 



in Eq. (||), but with the assumption of long-range cor- 
related disorder so as to stay within the soft-scattering 
regime. Specifically, the on-site disorder potential was 
defined by 



W X ^ 



where £ = ^ exp(— 2|r^|/£ 2 ), and fj were independent 
random variables drawn from a uniform distribution on 
[—1/2, 1/2]. The results are shown in Fig. El. 



• At zero energy the DoS vanishes in all cases but the 
extreme one of scatterers at the unitarity limit. In 
that particular case, spectral weight accumulates 
at the band center, reflecting the creation of im- 
purity bound states. For the special limit of zero 
chemical potential, corresponding to fully realized 
particle-hole nesting symmetry, the low-energy DoS 
diverges logarithmically in accord with the analysis 
of Lee and Pepin |l7j . 

• Away from zero energy (and for generic scatter- 
ers), a regime of linearly increasing DoS, tenta- 
tively identified as the DoS profile predicted by 
Senthil and Fisher fl3|| , is observed. We must cau- 
tion, however, that this identification does not con- 
vince us for weakly disordered systems: as dis- 
cussed earlier, the linear suppression appears in an 
insulating phase of separated localization volumes. 
Given the large size of the localization length in 
weakly disordered two-dimensional systems, it is 
not clear whether the separation of characteristic 
length scales can be realized on lattice sizes acces- 
sible to numerical computation. 

• Self consistency leads to a further suppression of 
the DoS, in particular in systems with binary- alloy 
type scatterers close to the unitarity limit. 

• As was explained in Ref . |86| , the nesting symmetry 
is of little relevance except for the case of unitary 
scatterers. 



B. Soft Scattering 

The quasi-particle spectrum of time-reversal invariant 
d-wave superconductors with soft scatterers (class All!) 
has been investigated numerically in Ref. [g5[. We will 
now review the results of that work in some detail. The 
starting point was the usual lattice Hamiltonian defined 
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FIG. 4. Density of states for A = 1, correlation length 
£ = 2, and disorder strengths W = 0, 1, . . . , 8, 10 (bottom to 
top at E = 0). All energies are measured in units of the hop- 
ping matrix element. The system size is L = 33 and the level 
broadening (introduced so as to suppress oscillations on the 
scale of the mean level spacing) F = 0.05. The inset shows 
the same data on a smaller scale with F = 0.0005. The finite 
DoS at E = is due to the finite level broadening F. 



The suppression of inter-node scattering by smoothing 
the potential is very strong: for a potential correlated 
over just two lattice spacings (£ = 2), the inter-node scat- 
tering matrix elements are reduced by a factor of about 
10~ 8 as compared to the intra-node matrix elements. Al- 
though such matrix elements will become relevant at very 
large scales, it is expected that potentials with £ > 2 
place small systems of size up to 100 x 100 firmly inside 
the pure symmetry class XIII. In practice, this means 
that each of the low-energy sectors associated with the 
four nodes in the dispersion relation is described by Dirac 
fermions subject to pure gauge disorder. 
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FIG. 
Fig. g 



5. Double logarithmic plot of the density of states of 
Disorder ranges from W = 1 to 10. Dots (•) repre- 
sent data and lines power law fits in the respective intervals. 
Inset: Density of states for W — 2 and L = 15 (dotted), 25 
(short-dashed), 35 (long-dashed), and 45 (solid). Note that 
the numerical uncertainties are considerably smaller than the 
amplitude of the fluctuations. 
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FIG. 6. Exponents a extracted from the fitted curves in 
Fig. U as a function of the disorder strength W for A = 1. 
The solid curve is the result of NTW, Eq. (ph. 



Figure || shows the DoS of the system for various disor- 
der strengths. The quantitative analysis of the spectral 
data shows that three different regimes can be distin- 
guished: (i) for low energies, e < -E m in, the structure of 
the spectrum is dominated by finite size effects. (For the 
lattice analysed in Ref. [^5| E min ~ E /10 where E de- 
notes the total width of spectrum.) The most apparent 
of these effects is a disorder and system size dependent 
bump in the DoS. At very low energies, the DoS van- 
ishes, as is seen in the high resolution inset of Fig. ||. 
(ii) For high energies, e > E max » E n /2, the struc- 
ture of the spectrum depends on non-universal lattice 
effects, (iii) Most interesting is the intermediate regime, 
E m - m < e < E max . In this region, the energy-dependent 
DoS exhibits power law behaviour (c.f. Fig. ||.) Of course, 
an energy window of width E m&x / E min ~ 5 provides a 
rather poor statistical basis for establishing power law 
behaviour. Nevertheless, the procedure seems justified 



as it is not just one power law with a single exponent but 
rather a two-parameter family of exponents a(W, 7) that 
is analysed. Here W measures the strength of the disor- 
der while 7 = t/A is the anistropy parameter. In Fig. [(], 
the exponents obtained by fitting the DoS determined 
numerically are compared with the NTW prediction ( p| ) 
for the isotropic case 7 = 1. Specifically, for the present 
system, 



p(E) ~ \E\ a , 



a 



I - 2g/ir 
l + 2.g/7r' 



9 = 
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Similarly, Fig. Q displays exponents obtained for fixed dis- 
order strength but different anisotropy parameters. No- 
tice that the comparison between the numerical data and 
the family of analytical exponents does not involve unde- 
termined fit parameters. The applicability of the NTW 
scaling law is limited to small disorder strengths, g < 1. 
In Ref. flsffl it has been argued that for larger values of g, 
the DoS becomes energy- independent. This prediction is 
supported by the numerical data. 

To summarize, numerical analysis of the quasi-particlc 
spectra in d-wave superconductors reveals the need to 
distinguish between three different types of disorder: 
scatterers at the unitarity limit, non- unitary point- like 
impurities, and soft scattering potentials. Although the 
comparison with analytical predictions is impeded by fi- 
nite size effects, there exists reasonable agreement for 
each of these types. 




0.01 



0.4 0.5 0.6 0.7 0.80.9 1 1.2 
E/A 



FIG. 7. Density of states for values of the order parameter 
A = 0.1 to 1.0 (top to bottom) and disorder strength W = 3. 
Each curve is shifted by a factor of 1.2 for clarity. 
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FIG. 8. Exponents a for W = 3 as a function of the order 
parameter A. The solid curve is the result of NTW, Eq. (fj3|). 



belonging to each nodal sector are locked by strong dis- 
order, the different WZW terms add and cancel pairwise. 

Finally, the formalism developed and investigated 
above is not entirely specific to d-wave superconductors. 
The global structures of the theory relied only on the ex- 
istence of a Dirac-like spectrum of the clean system. We 
believe that the general scheme outlined above could be 
applied in the investigation of other model systems with 
a gapless linear density of states such as gapless semicon- 
ductors, and superfluids. 
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VII. DISCUSSION 

This concludes our survey of the influence of disorder 
on the quasi-particle properties of disordered d-wave su- 
perconductors. Since we included the majority of the 
discussion in the introduction, we will limit our remarks 
to some key points: broadly speaking, the analysis above 
emphasized that the low-energy transport properties of 
the model d-wave system depend sensitively on the na- 
ture of the impurity potential. In two dimensions, a 
potential which is short-ranged in space places the sys- 
tem in the spin insulator phase, where all quasi-particle 
states are localized. On the other hand, a potential which 
contains only forward-scattering components leads to a 
marginally perturbed WZW theory in which the zero- 
energy quasi-particle states are critical, and the density 
of states vanishes as a power law. The low-energy theory 
in this case belongs to a one-parameter family of fixed 
points, each identified through a different value of the 
disorder-coupling strength. This is reflected in a low- 
energy density of states which varies with energy as a 
power law with an exponent that depends on the strength 
of disorder. 

Experimentally, the relevant scattering phenomenol- 
ogy is likely to be somewhere in between the cases con- 
sidered above: at intermediate energy scales, signatures 
of the critical theory may well be visible in thermal trans- 
port measurements, although the behaviour at very low 
energies must, ultimately, be that of the spin insulator. 

The collapse of the critical theory in the presence of 
strong disorder was attributed to a cancellation of the 
multi- valued WZW terms arising from the different nodal 
sectors of the theory. In fact, such cancellations are guar- 
anteed by symmetry to occur in lattice models which ex- 
hibit a Dirac-node structure at the Fermi level (such as 
the random 7r-flux model). In such lattice models, the 
nodes arise in pairs related by parity. Under the same 
parity transformation, the WZW term is mapped onto 
a partner with opposite sign. Therefore, when the fields 



APPENDIX A: GRADIENT EXPANSION AND 
THE CHIRAL ANOMALY 

To complement our derivation of the effective action for 
the soft-scattering limit from non-Abelian bosonization, 
we include here a derivation of the same action from the 
gradient expansion. The motivation is largely of ped- 
agogical nature: to facilitate comparison with existing 
works in the literature on weakly disordered fermion sys- 
tems, it is useful to present more than a single route 
to the construction of the critical theory. We also wish 
to point at some unexpected difficulties that are encoun- 
tered in the standard approach to deriving the low-energy 
effective action of the d-wave superconductor (or for that 
matter any disordered relativistic fermion system). 

To be specific we will formulate the gradient expansion 
for a soft-scattering system that is time-reversal invariant 
(class j4III). The inclusion of perturbations driving the 
model to any of the other three classes is straightforward. 
However, since this appendix mainly serves a pedagogical 
purpose, we will limit ourselves to the discussion of just 
one class. As discussed in Section EL the independence 
of the nodes entails a decoupling of the low-energy the- 
ory into two pairs of nodal sectors (1, 1) and (2, 2). We 
discuss one specific sector, say (1,1), anticipating that 
the full theory can later be obtained by straightforward 
combination of both sub-sectors. 

Our starting point is the soft-mode action for the (1, 1) 
sector derived in Section III and given by Eq. (p2]). Re- 



arranging matrix blocks, the action can be brought into 
the simpler form 



S[M] = STrln 



kM_~ 
d 



d 
kM 



(Al) 



where the block decomposition is in PH-space and we 
have omitted the superscript (1) on the derivative op- 
erator for notational simplicity. To further simplify the 
notation, we have set the two characteristic velocities Vi 
(i = 1,2) temporarily to unity (i.e. v = 1, and 7 = 1). 
(In fact, some authors attempt to get rid of these scales 
altogether by means of a coordinate rescaling xi — > ViXi. 
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However, as we are going to discuss below, this seemingly 
innocuous manipulation may lead to inconsistencies once 
the nodes are coupled. Moreover, the influence of such 
a rescaling on the unspecified source components of the 
action must be treated with caution. We will therefore 
re- instate the scales Vi towards the end of this section.) 
The bold-face notation STr means that we are taking the 
(super)trace over both superspace and Hilbert space. 

To compute a low-energy action from the above expres- 
sion, we can follow one of at least three different routes: 

• The most direct approach would be to introduce 
coordinates on the field manifold, say by M = c x , 
to expand around unity: M = 1 + X + . . ., and 
then to derive a low-energy action for the X's via a 
straightforward gradient expansion. Owing to the 
overall GL(2|2) invariance of the model, such an ap- 
proach determines the low-energy action not just in 
the vicinity of unity but rather on the entire man- 
ifold. 

• Alternatively, as in the main body of the text, one 
may resort to an entirely symmetry-oriented ap- 
proach and obtain the structure of the low-energy 
action by means of current algebra and non-Abelian 
bosonization. This was the route taken by NTW 

I- 

• Finally, a third option is not to introduce coordi- 
nates on the field manifold but to attempt a gradi- 
ent expansion directly based on the original degrees 
of freedom M. While such schemes are standard 
in applications of non-linear sigma models to disor- 
dered metallic systems, we here run into difficulties, 
caused by the appearance of ill-defined momentum 
integrals. The way to overcome this problem is first 
to subject (Al) to a UV-regularization scheme and 
only then to expand in the spatial fluctuations of 
the fields. This will be our method of choice in this 
section. Its main advantages are that it is com- 
putationally efficient and better exposes the global 
structures of the theory than a coordinate-based 
approach does. 



1. Chiral Anomaly 

Before subjecting the action functional to a gradient 
expansion, let us first make some pedagogical remarks. 
For the time being, let M be a field taking values in 
some matrix group G, and let us consider the functional 
determinant 



Z\M\ = 
D M = 



Bet D M 

kM' 1 
d 



TV In Da 



= e 

d 

kM 



-S[M] 



(When G is a group of supermatrices, Det has to be re- 
placed by SDet -1 .) Our goal is to expand ln(l/Z) in 
gradients to produce a low-energy effective action for M. 



Now, if G is a group of unitary matrices M 1 = , 
the determinant is not real: 



Z = Det 



kM 



-it 



-d 



so In Z has an imaginary part. On general field-theoretic 
grounds, we expect this imaginary part to be a multi- 
valued functional of WZW type. How can we compute 
Im In Z? 

A natural idea is to "take the square root": M = 
T 2 T^ , and manipulate the determinant as follows: 



Det( B nT 2Tl -^ = 



Det 



Ti 

T 2 ) \ T 2 ~ l dT 2 k J \ Tf 



One might now be tempted to assume the multiplicativ- 
ity of Det, which would lead to Z being equal to 



Z' = Det 



T 2 Y ~dT 2 



Dropping the diagonal factors on the left and right seems 
especially innocuous when T\ and T 2 are unitary. 

The motivation for trying to pass from Z to Z' is that 
the latter can be computed exactly by a standard proce- 
dure (see, e.g., Ref. Q and the next subsection) in the 
limit of small k. The result, 

lim lnZ'[Ti,T 2 ] = +2W[T 1 T 2 1 } , 
is expressed by the celebrated WZW functional: 

W[M] = — [ d 2 r Tr d u M~ x d u M + , 
1 1 16tt J p M 24tt ' 

T[M] = J d 3 r e^xTr M^d^M M~ l d v M M~ l d\M . 

Note, however, the inequality ReVF[M] > for uni- 
tary M. Thus, if Z were equal to Z' , the constant 
fields M (r) = Mq would minimize rather than max- 
imize the Boltzmann weight Z[M]. We would then 
be forced to conclude that the field theory with action 
S[M] — —lnZ[M] is unstable with respect to spatial 
fluctuations and does not exist. By extending the argu- 
ment to the supersymmetric setting, we would find the 
theory with action (Al) to be sick. On the other hand, 
we know (e.g. from non-Abelian bosonization) that this 
is not the case, so there must be something wrong with 
the present argument. Where is the error? 

The answer is that the manipulation taking Z into Z' 
disregards the existence of the notorious chiral anomaly 
and is correct only for T\ = T 2 , the case of a pure gauge 
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transformation. In other words, for a gauge transfor- 
mation with an axial component (T\ T2) the passage 
from Z to Z' is accompanied by a Jacobian different from 
unity. Indeed, for G = U(l), M = e llp , straightforward 
application of the method of Abelian bosonization Q 
gives 

-lnZ[e^] = ±- [ d 2 r(d^) 2 = lnZ'[e^] . 
oir J 

By analogy, we expect that also in the non- Abelian case, 
correct evaluation of Z[M] yields a stable theory with the 
proper sign of the coupling. A safe way of computing the 
gradient expansion is to first UV regularize the Dirac op- 
erator Dm after which axial gauge transformations can 
readily be performed. 

Returning to our original problem, we notice that a 
technically convenient way o f re gularizing in the ultravi- 
olet is to add to the action (Al) a term 



STr In {^ M B 





eM 



which vanishes by supersymmetry when e is taken to be 
a positive infinitesimal (e — * 0+). The resulting expres- 



the same procedure was applied to the non-Abelian 
Schwinger model, i.e. massless 1 + 1 dimensional Dirac 
fermions coupled to an SU(iV c ) gauge field.) Without 
loss, we equate T = T\ = T% and rewrite the action in 
the form 



sion, 



( e T~ 1 dT 

= STr In (e 2 - TdT^&T) , 

where a UV cutoff at the momentum scale k is implied. 

The proper-time regularization scheme for an elliptic 
operator H is implemented by 

r°° 

-TrlniT= / — Tre"^ . 

Notice that the lower integration bound 1/fi cuts off the 
contributions to InDetiJ from eigenvalues of H greater 
than fi and thus regularizes in the ultraviolet. Applying 
this scheme to our action (with cutoff /j = k 2 ), we obtain 

S = r - STr e s(e 2 -TBT-^T) 

Jl/K* S 



S[M] = STr In 



kM- 1 

8 



d 

kM 



zM_~ 
d 



B 

zM 



is indeed manifestly well-behaved in the ultraviolet. (The 
difference between eM and kM becomes negligible for 
large eigenvalues of the Dirac operator, in which case the 
two matrix factors cancel each other and the action ap- 
proaches zero. ) Setting M = T 2 T~ Y and using the cyclic 
invariance of the trace, we rewrite the action functional 
as 



STr In 



STr In 



Tx X dT x 



e_ Tf 1 ^ 



T^dT 2 



Because k now acts as a mass, the low-energy limit of 
the theory is captured entirely by the second term. The 
first contribution becomes appreciable only for momenta 
larger than k, where it cancels the second term. Thus 
the role of the first contribution has been relegated to 
that of a UV regulator. We are, of course, at liberty to 
replace it by some other UV regularization scheme. Do- 
ing so and expanding the second term in gradients, we 
safely arrive at the WZW action, now with the correct 
overall sign. This will be demonstrated in more detail in 
the next subsection. 



2. Heat kernel regularization 

For our purposes, it is convenient to use Schwinger's 
proper-time regularization (see, e.g., Ref. Q where 



This expression is both UV and IR finite and could in 
principle serve as the starting point for a gradient ex- 
pansion. Much easier than the direct evaluation of S, 
however, is the evaluation of its variation SS. We will 
therefore proceed by varying S with respect to some pa- 
rameter, say t; then we will compute SS = S, and finally 
we will reconstruct S by integrating S with respect to t. 

Thus we consider some one-parameter family of fields 
T(r,t) with T(r,Q) = 1 and T(r, 1) = T(r), and we dif- 
ferentiate with respect to t. This results in 

S= f ds STr (fdT- 2 dT + TdT- 2 df 

- Td(T- lr tT- 2 + T- 2 TT~ 1 )dTje~ a ( s2 - T8T ~ 2dT \ 

We next use the cyclic invariance of the supertrace to 
convert the integrand into a total derivative with respect 
to the integration variable s: 

/>oo 

S= ds e~ s2s STr (T _1 T + TT _1 ) 

Jl/K 2 

x j ^ e sTdT' 2 dT _ e sT~ 1 dT 2 dT^ 1 ^j 

Performing the integral over s, setting the infinitesimal 
e to zero, and making the integration over real space 
(J d 2 r) explicit, we obtain the expression 

S= [ d 2 r STr (T^f + TT- 1 )(r) 



x . j, c -K H-L _ e -K H 2 r \ 
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where 

Hi = -rtofSoT 1 , 
H 2 = -TdoT- 2 doT , 

and the symbol o means composition of operators. 

Our next task is to compute the diagonal parts of the 
heat kernels (Y\e~ sHi \r') (i = 1,2), for small values of 
the dimensionful parameter s — 1/k 2 . This is a standard 
exercise in semiclassical analysis, and its solution can be 
found in textbooks BG]. Re-expressing Hi in the form 



Hi 



4(0„ - iA.f + B 



where the non-Abelian gauge potential and field 
strength B are functions of T and its derivatives (which 
for brevity we do not specify here) , we have the standard 
short-time expansion 



1 

7TS 



B(v) 



0(s) 



The same can be done for H 2 instead of Hi. By taking 
the difference of the two expansions, we obtain 



[T _1 dT, dTT- 1 ] 



+d(dTT- L ) + d(T- L dT)\ (r) + 0(1/ kT) . 

On dimensional grounds, the term 0(1/ n 2 ) must involve 
four derivatives and therefore becomes negligible for wave 
lengths much larger than the short-distance cutoff 1/k. 

The above expansion is now substituted into the ex- 
pression for S. We then arrive at 

S = — J d 2 rSTr(T- 1 f + TT- 1 )(d(T- 1 dT) 

+ d(dTT- 1 ) + T~ 1 dTdTT~ 1 - dTT~ 2 dT 

It is not hard to verify that, on making the identification 
M = T~ 2 , this expression coincides with 

1 d f 

S = / (frSTriduMduM- 1 ) 

+ i^J ^e^STr^tMAf-^A/Af-^MAf- 1 ) . 

Integrating over time and noticing that S = f Q Sdt, we 
obtain the WZW functional given in (|30|), with g = 0. Fi- 
nally, we undo the rescaling made at the beginning of the 
calculation, and arrive at the anisotropic effective action 

S[M, 7 ] = -UiM]- J- j d 2 rx 

x STr (j^diM^diM + 1 - 1 d 2 M^ 1 d 2 M) . 

The corresponding result for the other pair of nodes (2, 2) 
is obtained by exchanging coordinates Xi <-> x 2 . 



Notice that this result has the peculiar feature of being 
independent of the disorder. In the isotropic case 7=1, 
the model becomes completely universal in the sense that 
its coupling constants assume fixe d values. By the argu- 
ments reviewed in Appendix B 2 , this WZW model is 



equivalent to free relativistic fermions, i.e. our original 
model in the absence of disorder. But this raises the 
question where the information on the presence of impu- 
rities was lost. After all it is hard to conceive that the 
disorder strength g should be reflected only in the value 
of the UV momentum cutoff k ~ e~ 47rt A ^ 9 . 

Our analysis using non-Abelian bosonization described 
in the main text reveals the fact that the action above 
should be supplemented by a current-current interac- 
tion with coupling given by g. This begs the question 
how the existence of this term got lost in the standard 
scheme of saddle-point approximation plus gradient ex- 
pansion. The key to the answer of this question lies in 
the presence of massive modes in the model, a fact we 
have ignored thus far. Indeed, we had immediately re- 
duced the Hubbard-Stratonovich field Q to its Goldstone- 
mode content, Q — > inM . In fact, however, the field Q 
also contains massive modes, i.e. modes P which are not 
compatible with the chiral symmetry of the Hamiltonian 
and, therefore, fluctuate at a finite energy cost. One way 
of handling the situation would be to put Q — inPM , 
where the P's are not just set to unity but integrated out. 
This procedure results in the appearance of an additional 
current-current interaction 

STr(Af- 1 aM) STi(MdM^ 1 ) . 

Unfortunately, it turns out to be impossible to reliably 
determine the value of the coupling constant of this per- 
turbation within the standard scheme: to obtain the ad- 
ditional term, one has to integrate out massive fluctu- 
ations, and yet the mass gap characterizing the modes 
P does not suffice to justify a Gaussian approximation 
to this integral. The reason for the last fact is that 
the disorder-generated perturbation of the field theory is 
strictly RG marginal, which means that there is no dy- 
namically generated mass scale in the problem. Thus, no 
intrinsic mechanism stabilizing the Gaussian approxima- 
tion exists, fluctuations are important, and to determine 
the coupling constant, the P-integral must be performed 
exactly. It goes without saying that this is difficult to do 
in practice. Ultimately it is this deficiency of the stan- 
dard scheme which forces us to build our theory on the 
less standard approach of non-Abelian bosonization. 



APPENDIX B: DIRAC FERMIONS IN A 
RANDOM VECTOR POTENTIAL 



As was reviewed in the main text, the low-energy 
quasi-particles of a dirty d-wave superconductor behave, 
in the single-node approximation, as Dirac fermions in 



3G 



a random vector potential. We have argued that non- 
Abelian bosonization takes that theory into a supersym- 
metric WZW model, from which we rederived the crit- 
ical exponent for the density of states. In the present 
appendix we elaborate on this issue and supply more of 
the technical details. 



1. WZW model of type A\A 

Let us begin by establishing the general context with 
a few remarks. The characteristic feature of a WZW 
model is the multi- valued term T[M]. Multi- valued func- 
tional of this type were studied in the context of Hamil- 
tonian mechanics by Novikov pl[ , and became firmly es- 
tablished in field theory through Witten's celebrated pa- 
per |3f| on non-Abelian bosonization. Ever since, WZW 
models defined on compact groups have been an object 
of intense study. There exists a vast amount of litera- 
ture on them, and they have been solved in great de- 
tail. To a much lesser extent, field theorists have also 
studied WZW models of the non-compact type. The tar- 
get spaces of these models are not groups but are non- 
compact symmetric spaces, the simplest example being 
SL(2, C)/SU(2). One of the rare occurrences of such a 
model is found in Q . 

The target space of the WZW model we will be con- 
cerned with transcends the classical setting in that it 
is a superspace. It turns out that the proper mathe- 
matical construction of a WZW model with superspace 
target is a delicate matter. Indeed, to define a func- 
tional integral on Euclidean space-time, one needs a tar- 
get space with a Riemannian structure, providing for an 
action functional that is bounded from below. Sadly, the 
invariant geometry on supergroups, such as GLc(n\ri), 
U(n|n), GL n(n\n), OSp R (2n\2n) etc., or even on sym- 
metric quotients such as GLc(n\n)/XJ(n\ri), is never Rie- 
mannian, but always of indefinite signature. Therefore, 
WZW models, and non-linear sigma models in general, 
do not exist on supergroups, at least not in the literal 
sense (i.e. without some procedure of analytic continua- 
tion of the fields). 

One can easily appreciate this point by looking, for 
example, at the complex Lie supergroup GLc(l|l), 
with the standard (bi-)invariant metric given by k — 
-STrdg-Mff, g G GL C (1|1). To understand the prop- 
erties of this metric tensor, it suffices to examine it on 
the tangent space at the group unit, which is the com- 
plex Lie algebra gl c (l|l). Elements of this Lie algebra 
are written as 

with commuting a, b and anticommuting a, f3. Via the 
exponential mapping g = e x , the invariant metric k re- 
stricts to 

k\ = STr dXdX = da 2 - db 2 + 2dadj3 . 



Clearly, the numerical part da 2 — db 2 is of indefinite sign, 
and this does not change on passing to gl fi (l|l) by re- 
stricting the complex number field to the reals. The same 
remains true if we choose a compact real form u(l|l) by 
setting a — itpi and b — iip2, with real yi,2- Thus no mat- 
ter which real form of gl c (l|l) we choose, the numerical 
part da 2 — d& 2 of the invariant metric will always be of 
indefinite sign. Note also that it is pointless to replace 
the compact supergroup U(l|l) by its non-compact ana- 
log GLc(l|l)/U(l|l). This just reverses the overall sign 
of the metric, leaving it as indefinite as before. 

One might object that if the target space were taken to 
be the compact supergroup U(l|l), the range of all field 
variables would be finite, and one could then say that 
compact integrals always exist, no matter what is the sign 
of the exponent of the Boltzmann weight. This reasoning 
would be valid if we were dealing with a finite number 
of compact integrals. However, the system at hand is a 
functional integral over infinitely many degrees of free- 
dom, and to ensure the existence of this functional inte- 
gral, we do need a target space with Riemannian struc- 
ture. Indeed, if the metric of the Euclidean field theory 
were Lorentzian rather than Riemannian, the energy of a 
field configuration could be diminished below any bound 
by shortening the wave length of some field. For the case 
of U(l 1 1) with Boltzmann weight exp — j STr(g^ 1 d fJ ,g) 2 
this would be the field a = iipi . If one attempted to define 
the field theory by its Feynman diagrams, one would en- 
counter the problem that the free-field approximation to 
the Boltzmann weight for a fluctuation <pi(k) with wave 
number k is exp(|fc| 2 |</?i(/c)| 2 ), which increases with the 
strength of the fluctuation! Thus the field theory would 
be unstable w.r.t. short-wave length fluctuations in ipi. 
Of course one could rescue the stability by putting the 
theory on a lattice with, say, nearest-neighbour coupling 
for the U(l|l) field. However, the state of lowest energy 
of such a lattice theory would be ferromagnetic in if2, 
but antiferromagnetic in the ip\. As a result the ground 
state would not be invariant under global supersymme- 
try transformations. Even worse, the partition function 
of the theory (in the absence of sources) would be iden- 
tically zero and hence not normalizable to unity. This 
follows essentially from the fact |)3| that U(l|l) has van- 
ishing volume w.r.t. its invariant Berezin-Haar measure. 
All these considerations make it rather unlikely that one 
could ever make sense, for our purposes, of a Euclidean 
space-time nonlinear sigma model or WZW model with 
target space U(l|l) or more generally U(n|n). 

The construction that solves the difficulty was de- 
scribed in a general mathematical setting in Ref . |5(J . Its 
field-theoretic implementation was discussed in some de- 
tail in Section 6 of Ref. @ and Section 7.1 of Ref. g§. 
The basic idea is easily understood. To make the nu- 
merical part of the metric tensor da 2 — db 2 positive def- 
inite, we should take a = x from the real numbers, and 
b = iy from the imaginary numbers. By returning to 
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group level via exponentiation, we then get manifolds 
M F = U(l) and M B = GL(1, C)/U(l). In the gen- 
eral case of n Green functions, we obtain M F = U(n) 
and M B = GL(ro, C)/U(n). These are Riemannian sym- 
metric spaces of compact and non-compact type, respec- 
tively. On incorporating the fermions into the exponen- 
tial mapping, we are led to some sort of superspace. An 
equivalent procedure is to start from the complex su- 
pergroup GLc(n\n) and restrict the bosonic degrees of 
freedom to M B x M F . (This construction does not let 
us impose any reality constraints on the fermions [ p4| . 
Fortunately, there exists no fundamental principle that 
would force us to do so.) Let X„ denote the resulting 
space. X„ is not a group, but belongs to the category 
of Riemannian symmetric superspaces |50| . It is called 
type A\A, which tells us that both the BB and the ff 
sector are symmetric spaces belonging to the A (or "uni- 
tary") series. Its main property is that the metric tensor 
k = — STr dg~ 1 dg restricts to a Riemannian structure on 
M B x M F . 

Thus, the field theory to be studied has the target 
space X„. We denote the field by M, and set n = 2 
(although nothing essential depends on that choice). By 
slight generalization of the action functional ( |30| ) , we con- 
sider 

S[M] = kW[M] 
h\ r 

+ — I d 2 rSTr(MaM~ 1 )STr(Af- 1 aM) 

with W[M] being the WZW functional for GL(2|2): 

W[M] =--!-/ (frSTrd ll M- 1 d IM M+^^ . 

This is the theory we shall now analyse. As a by-product 
we will get (by setting A = 0) a justification of the 
bosonization rules that lead to (|3C|). 

2. Functional integral solution 



The main tool for dealing with the WZW functional is 
a relation due to Polyakov and Wiegmann |)5| : 

W[gh] = W[g] + W[h] + - [ d 2 r STr g^dg dhh~ l . 



Using it, one sees that W[M] is invariant under local 
transformations 

M(z,z) i-> g(z)M(z,z)h(z)~ 1 , 

where g(z) and h(z) take values in GLc(2|2). We will 
see that a deformation of this symmetry survives in the 
presence of a no n- vanishing coupling A. 

Consider now the expectation value {T[M\) where 
!F[M) is some functional of the field M concentrated 



at a set of points ri, T2, rjy. Exploiting the invari- 
ance of the functional integration measure under left 
translations, we make a change of integration variables 
M(r) h- > e _x ^ r ^M(r), and denote the first variation of 
the field by 6 x M(r) = -X(r)M(r). The resulting vari- 
ation of the action defines the current J: 



S X S = -— I d 2 rJ 5x 



Jx = -kSTiiXMdJVr 1 ) 



2k\ 



STr(X)STr(MdM- 1 ) 



To derive this expression for the current, we made use of 
the Polyakov- Wiegmann formula, and we used the iden- 
tity 

aSTr(M~ 1 9Af) = BSTriM^dM) 

in conjunction with two partial integrations to cancel 
some terms. (The boundary terms produced by these 
partial integrations at the points r±, . . . ,rjy cancel each 
other.) 

The invariance of the expectation value (J-[Mj) w.r.t. 
the variation S X M implies 

(6 X T[M\ + F[M]^J d 2 r J dx ) = . 

From this relation it immediately follows that 8Ja — 0, 
called the equation of motion, holds under the functional 
integral sign (away from the points of support of JF[il/]) 
for any spatially constant A £ gl c (2|2). As an additional 
consequence one has 



1 

2vri 



dz(J x (z)-F[M}) = {8 X T[M}) , 



where 7 is any closed contour that circles once around the 
points of support of JF[M]. By specializing to ,F[M] = 
M(0)!Fi[M], where T\ is some other functional, one gets 
the operator product expansion between the current and 
the field: 



Ja(z)M{0) = - 



AM(0) 



and by taking F[M] = J B (0)J c 'i[M] one gets the OPE 
for the currents themselves: 



Ja{z)Jb(0) 



kf(A,B) , J[A,B](°) 



2A, 



f(A, B) = -STr(AS) + — STr {A) STr (B) . 

The dots indicate terms that remain finite in the limit 
z -> 0. 

Identical considerations can be made starting from 
field variations <VM(r) = M(r)F(r), which generate 
right translations M(r) 1— ► M(r)e y ^ r ^. They lead to a 
conserved current 

J Y = —kSTr(YM~ 1 dM) + STr(F)STr(M- 1 9M) , 

7T 
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which is antiholomorphic, i.e. satisfies 8Ja = for spa- 
tially constant A. 

The behaviour of correlation functions under confor- 
mal transformations is determined by the stress-energy- 
momentum tensor, written in components as 



Taydx^dx v = T zz dz 2 



Conformal invariance implies T zz = = T zz , and dT zz = 
= dT zz (on solutions of the equations of motion) . We 
focus on the holomorphic sector and set T(z) = 2irT zz (z). 
Classical considerations based on the general formula 
T^v = (df,(p l )dC/d(d u <f l ) - S^C suggest 

k k\ 
Tci = -STr (MdM^ 1 ) 2 STt (MdM ) , 

2 7T 

but this, it turns out, is modified by quantum fluctu- 
ations. The correct expression for T(z) is obtained by 
demanding that the leading singularity in the operator 
product T(z)Ja(0) be z~ 2 Ja(0), as results from J a be- 
ing a holomorphic conserved current. 

The explicit form of the stress-energy-momentum ten- 
sor is obtained by the Sugawara construction, which rep- 
resents T(z) as a quadratic form in the currents. Let {e,} 
be some basis of the Lie superalgebra gl(2, 2) with metric 
Kij = — STre^e.,. Indices are raised by k ij Kju = 5\. We 
write Ji = J Ci for short, and J e = J;d for the current 
corresponding to the unit matrix. Setting 

f(z) = ±-K i >:J i (z)J j (z): 

where the colons mean normal ordering (i.e. subtraction 
of the short-distance singularities of the operator product 
at coinciding points), and using the OPE for the currents 
and the associativity of the operator product algebra, we 
get 



lim z 2 T(z) J A (0) 



2A 

= J A (0)-— STr(A)J e (0) 

7T 

+ ^T J [e*,[e,',A]](0) ■ 



The last term is a quantum correction, which vanishes 
in the limit fc — > oo. For gl(2,2) one verifies the rela- 
tion /e* 3 '[ei, [ej, A]] = 2eSTrA This allows to rewrite the 
right-hand side of the preceding relation as 



Ja(0) 



1 2A s 

_ ] STr(^)J e (0) 



To obtain the required form of the operator product 
T(z)Ja(0) — Ja(0)/z 2 + we must cancel the second 
term, which is achieved by defining 

T(z) = f(z) + 1 ~ 2 2 fc A 2 fc/?r : Je(z)J e (z) : . 

On setting k — 1 and renaming A to g, we arrive at the 
expression for T(z) claimed in Eq. (p3). Note that the 



symbol for normal ordering can be dropped, as the ex- 
pression is already finite as it stands |39[] . 

Our final step is the computation of the OPE of T(z) 
with the fundamental field M(0). This is readily done 
from the formula for T{z) and the OPE between the 
currents and M(0). Using the fact that the quadratic 
Casimir invariant evaluated in the fundamental repre- 
sentation of gl(2,2) vanishes (n^eiej = 0), we find 



■yZ 



where Am = (1 — 2Afc / V) / '(2fc 2 ) is the (holomorphic) 
scaling dimension of the field M. The total dimension is 



1 - 2Afc/7T 



Note that the dimension becomes negative for disor- 
der strengths A greater than n/2k. If we substitute 
a a 



M = exp 



/3 b 



into the action S[M] and isolate the 



terms involving the boson-boson field a = X BB only, we 
get 



S[M] = 



1 - 



2A 



We see that the coupling constant turns negative, so the 
functional integral becomes unstable, at A = 7r/2. For 
k = 1 the onset of instability is located right at the value 
of A where the dimension of M becomes negative. A 
plausible scenario for what happens beyond that point 
(for k = 1) has been suggested by Gurarie |p7| . For 
higher level k > 1, there exists an intermediate regime 
n/2k < A < 7r/2, where Am + Am is negative while the 
functional integral is still stable. In this regime we ex- 
pect the density of states (the expectation value of the 
field M) to diverge at zero energy. 

For k = 1 and A = 0, the conformal dimensions of 
M are (Am, Am) = (1/2, 1/2). These coincide with the 
conformal dimension of the bilinears ^i^j in the super- 
symmetric Dirac theory with Lagrangian 

C = + *f<9*i . 

Since the OPEs involving the currents match, too, we ex- 
pect that the two theories are equivalent. The bosoniza- 
tion rules are the same as in (g4|) but for factors of two: 

*J*f <-> TT^Mdhr 1 , 



dM 
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APPENDIX C: WZW THEORY FOR SYSTEMS 
WITH BOUNDARY 



In the main text we had said that a single WZW the- 
ory represents too narrow a framework to consistently 
describe systems with boundaries: the WZW functional 
exp iT[M]/12ir for maps M from the position space E 
to the target space Q (the Riemannian symmetric super- 
space of type A\A in GL(2|2)) makes sense only if E is 
closed (i.e. has no boundaries.) The reason is that, to 
define T[M] as an integral of the 3-form STr(M _1 dM) A3 , 
one needs to extend M to a map M : B — > Q where B is 
some 3-space. But in order for the WZW functional to be 
independent of the choice of extension M, the boundary 
of B must coincide with E (dB — E), which is possi- 
ble only if E is closed (<9E = ddB = 0). Thus the WZW 
functional is ill-defined in the presence of a boundary. To 
the extent that the WZW model is a faithful representa- 
tion of the Dirac theory, we expect the implementation 
of a boundary to be problematic in the latter, too. 

How can we cure this disease and arrive at a mean- 
ingful functional integral for systems with a boundary? 
There exists only one solution to the problem f8l"ll , and it 
requires having a minimum of two WZW theories. Let's 
assume there exists a pair of WZW fields M and M', with 
the sum of WZW terms being T[M] +T[M'}. The bound- 
ary condition we impose is M _1 j as = M'| 0E . (To ensure 
current conservation, we also need to impose a condition 
on the normal derivatives of the fields at the boundary.) 
This boundary condition guarantees that the two maps 
M : E — > G and M' : E — > Q combine to a single map 
E U (— E) — * Q, where E U (— E) consists of two copies of 
the position space, glued together at the boundary <9E. 
Because the double covering EU(— E)isa closed man- 
ifold, we can regard it as the boundary of a 3-space B, 
so the independence of the choice of extension is restored 
and the total WZW functional is again well-defined. 

Having understood this, the question arises: where do 
we take the second WZW theory from? Depending on 
the physics of the problem at hand, a number of options 
are conceivable: If the original problem mapped onto an 
even number of WZW species anyway - this is the case 
for our d-wave system - it is natural to glue these fields 
pairwise together at the boundaries. This procedure is 
discussed in the text. However, one may also contem- 
plate situations where the low-energy sector of the orig- 
inal problem contains only a single WZW theory. (In 
the language of fcrmions, this amounts to considering a 
single Dirac fermion species.) In this case, invocing a sec- 
ond WZW field or, equivalently, a second Dirac species, 
changes the low-energy content of the theory. To repair 
this, we must give a mass to the second species, so as to 
prevent it from contributing to the long-range behaviour 
of the correlation functions. Note the important fact that 
the mass of the second species breaks parity, so we can no 
longer conclude on symmetry grounds that pseudoscalar 
observables such as the Hall conductivity vanish! If wc 



again carry out the low-energy reduction to the nonlinear 
sigma model, we will get the same answer as in section 
IV D, as the auxiliary second species drops out in the 



low-energy limit. Moreover, we expect the sign ambigu- 
ity &xy = ±1/2 to be resolved: it is the sign of the mass 
of the second species that will determine the sign of a xy . 

The second scenario is not new but is known from the 
work of Ludwig et al j73j who studied a disordered elec- 
tron model on a square lattice with a magnetic flux of 7r 
per plaquette. The low-energy limit of that model con- 
sists of two species of Dirac fermions, one with a light 
mass and the other with a heavy mass. The Hall con- 
ductivity of the model does not vanish when the light 
species becomes massless (and hence the low-energy sec- 
tor parity- invariant), and is determined by the sign of the 
Dirac mass of the heavy species. 
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